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SUMMARY 


The  scattering  of  sunlight  by  interplanetary  dust  gives  rise  to  the 
zodiacal  light,  a  ubiquitous  feature  of  the  night  sky  which  is  the  limiting 
background  for  most  IR  observations.  This  changing  (with  wavelength,  look 
direction,  and  distance  from  the  sun)  background  can  be  understood  only  in 
the  larger  context  of  the  physics  of  interplanetary  dust.  We  have  obser- 
vationally  separated  the  zodiacal  light  from  other  astronomical  sources, 
devised  a  mathematical  inversion  to  extract  the  maximum  amount  of  informa¬ 
tion  about  the  dust  from  space  observations,  observationally  proved  that  the 
dust  complex  is  neither  homogeneous  nor  simply  distributed  through  the  solar 
system,  and  predicted  that  the  dust  mav  tend  to  accumulate  in  enormous  arcs 
which  span  the  solar  system.  The  inversion  technique  is  the  most  promising 

method  to  extract  from  visual  and  near  IR  observations  the  parameters  neces- 

i 

sary  to  model  the  IR’  emission  of  interplanetary  grains. 

OBJECTIVES 

The  scattering  of  sunlight  by  interplanetary  dust  is  referred  to  as 
zodiacal  light,  since  the  brightness  and,  therefore,  the  dust  distribution 
have  a  maximum  near  the  zodiacal  band.  Only  the  zodiacal  light  contains 
information  on  the  integrated  properties  of  the  whole  ensemble  of  interplane¬ 
tary  dust.  Observations  of  brightness,  polarization,  and  color  contain  in¬ 
formation  about  the  size,  shape,  composition,  and  overall  optical  properties 
of  the  grains.  Observations  made  from  different  locations  in  the  solar  sys¬ 
tem  provide  data  on  the  change  in  zodiacal  light  with  heliocentric  distance 
and,  therefore,  on  the  spatial  distribution  of  the  dust  and  on  changes  in  its 


1 


optical  properties  as  functions  of  sup-particle  distance  and  height  above  or 
below  the  ecliptic  plane. 

We  have  accumulated  a  vast  l.brary  of  unique  observational  data  on  the 
brightness,  polarization,  and  color  of  the  zodiacal  light,  including:  exten¬ 
sive  ground-based  observations  of  the  light  of  the  night  sky  and  of  Comer 
Ikey-Seki  (1965  VIII)  from  Mt .  Haleakala,  Hawaii,  the  first  multicolor  obser¬ 
vations  of  large  regions  cf  the  sky  from  space  (Skvlab) ,  and  the  first  obser¬ 
vations  of  large  regions  of  light  with  heliocentric  distance  (Pioneers  10  and 
11).  Each  of  these  data  sets  is  being  used  with  laboratory  and  theoretical 
studies  to  provide  information  on  the  spatial  distribution,  optical  and  physi¬ 
cal  properties,  and  dynamics  of  the  interplanetary  dust. 

In  October  1976,  we  began  an  Air  Force-funded  program  to  complete  the 
reduction  and  analysis  of  the  Pioneer  10/11  and  Skylab  data  and  to  combine 
these  data  with  ground-based  results  and  theoretical  studies  to  map  zodiacal 
light  brightness  and  polarization  over  the  sky  and  to  determine  the  nature  and 
distribution  of  the  interplanetary  grains.  Studies  have  been  made  of  the  dy¬ 
namics  of  small  particles  in  the  solar  system,  and  of  methods  for  inverting  the 
zodiacal  light  brightness  integral  to  derive  information  on  the  spatial  dis¬ 
tribution  and  optical/phvsical  properties  of  the  grains. 

The  fundamental  goal  of  our  research  —  to  define  the  spatial  distribu¬ 
tion,  the  sinks  and  sources,  the  dynamics,  and  the  physical  description  of 
the  interplanetary  grains  —  is  a  two-step  process.  The  first  step  consists 
of  observing  and  defining  the  macroscopic  features  of  the  dust  complex.  The 
second  step  (now  in  its  beginning  stages)  is  to  use  the  physics  of  light 
scattering  to  derive  the  microscopic  features  of  the  grains.  During  this 
fourth  year  of  AFOSR  sponsorship,  we  have  established  step  one  firmly  enough 
to  be  used  as  the  "‘"“pping  stone  to  step  two.  The  accomplishments  of  the 
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i;isl  year  in  capsii'nim:  sup  one  are  briefly  listed  below, 
re! erences  c 1  lea  t Herein  are  results  01  previous  ArOSK  sponsorship  and  are 
provided  tor  completeness.  Publications  fully  ui  partially  supported  by 
this  contract  are  included  as  attachments. 


SIGNIFICANT  RESULTS 


1.  Background  starlight  over  the  sky.  Pioneer  10/11  observations  from 

beyond  the  asteroid  belt,  where  the  zodiacal  light  was  found  to  be 

negligible,  provided  the  first  all-sky  data  on  background  starlight 

in  the  absence  of  the  atmospheric  radiations  and  zodiacal  light.  These 
two-color,  all-sky  data,  now  in  the  form  of  an  Atlas,  are  a  fundamental 
source  of  data  for  a  number  of  fields  in  astronomy  and  astrophysics. 
(Schuerman,  et  al.,  COSPAR  1977).  Sample  data  are  illustrated  in  a 
recent  issue  of  Sky  and  Telescope  (February  1981,  attached),  and  the 
full  body  of  data  is  being  considered  for  use  with  a  new  edition  of 
the  well-known  Atlas  Coeli. 

2.  Inversion  of  the  zodiacal  light  (ZL)  brightness  integral.  The  general 

inversion  of  the  ZL  brightness  integral  has  been  derived  and  is  being 

used  in  our  analysis  of  Pioneer  ZL  observations.  It  is  also  in  use  or 
planned  for  analysis  of  other  ZL  observations,  ground-based  or  Earth¬ 
orbiting  or  space  probe.  This  is  a  major  tool  and  the  only  complete 
source  of  information  on  particle  properties  from  photometric  data. 
(Schuerman,  Bull.  AAS,  _9,  564,  1977;  Schuerman,  Planetary  Space  Sci., 

27,  551,  1979;  Schuerman,  Planet.  Space  Sci.,  28,  1077,  1980.) 

3.  Properties  of  interplanetary  dust  beyond  1  AU.  Our  inversion  of  Pioneer 
10/11  ZL  data  has  found  the  first  evidence  that  the  properties  of  dust 
beyond  1  AU  are  not  homogeneous  —  one  of  the  major  assumptions  in  most 
studies  of  the  interplanetary  dust  complex.  (Schuerman,  in  Proc.  IAU 
Sump.  So.  90,  Halliday  and  McIntosh,  eds.,  D.  Reidel  Publ.  Co.,  Dordrecht, 
25 /  1980.) 

4.  The  positions  of  maximum  interplanetary  dust  density.  Our  analysis  of 
ground  and  space  observations  of  the  photometric  axis  of  zodiacal  light 
has  shown  that  there  is  no  unique  plane  of  "symmetry"  or  concentration 
for  interplanetary  dust,  but  rather  a  "multiplicity  of  planes"  apparently 
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associated  vith  the  orbital  planes  ot  the  planets.  (Misconi  and  UVinbeig, 
Z2ZCK?c,  200,  IAS-,,  1978;  Misconi,  in  Fvza.  ZAl'  .-'a.rr .  el,  49.  1980. 

5.  As ymm.e tries  in  the  spatial  distribution  of  interplanetary  dust.  As  part 
of  a  study  to  include  radiation  pressure  and  the  Poynt i ng-Robertson  effect 
in  the  classical  restricted  three-bodv  problem,  it  vis  discovered  that 
thexe  might  exist  interplanetary  dust  arcs.  This  leads  to  a  number  of 
ramifications  concerning  sources  and  sinks  for  the  interplanetary  dust 
complex.  (Schuerman,  Asiro’p'r'.uc .  J. ,  238 ,  337-342,  1980;  Schuerman,  in 
lO’x c .  IAU  Jxro.  i'C .  9G,  285,  1980. 

6.  Ph.D.  dissertations.  Three  theses  have  been  completed: 

(1)  G.  N .  Toller  —  A  Study  of  Galactic  Light,  Ext raga lac t i c  Light,  and 
Galactic  Structure  Using  Pioneer  10  Observations  of  Background  Star¬ 
light  (Space  Astronomy  Laboratory  and  SUMY  Stony  Brook). 

(2)  R.  Schaefer  —  Calculations  of  the  Light  Scattered  by  Randomly 
Oriented  Ensembles  of  Spheroids  of  Size  Comparable  to  the  Wavelength 
(Space  Astronomy  Laboratory  and  SUNY  Albany) . 

(3)  B.  S.  S.  Gustafson  —  Scattering  by  Ensembles  of  Small  Particles: 
Experiment,  Theory,  and  Application  (Space  Astronomy  Laboratory  and 
University  of  Lund). 

Their  abstracts  are  attached. 
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THE  RESTRICTED  THREE-BODY  F’ROBLEM  INCLUDING  RADIATION  PRESSURE 

Donald  W.  Schlerman 

Space  Asironomv  Laboraiorv.  State  University  of  New  York  at  Albany 
Reeened  /V**V  September  14  accepted  /V’%'  Ywt'm/vr  14 

ABSTRACT 

The  classical  restricted  three-body  problem  is  generalized  to  include  the  force  of  radiation 
pressure  and  the  Poynting-Robcrtson  effect.  The  positions  of  the  Lagrangian  points  F4  and  L<  are 
found  as  functions  of  ft.  the  ratio  of  radiation  to  gravitational  forces.  The  Poynting-Robcrtson 
elfect  renders  the  F4  and  F<  points  unstable  on  a  time  scale  ( T)  long  compared  to  the  period  ol 
rotation  of  the  two  massive  bodies.  For  the  solar  system.  Fis  given  by  Fs  [(1  -  ft)2  '  /f]544u:  yr. 
where  a  is  the  separation  between  the  Sun  and  the  planet  in  A  Li.  Implications  for  space  colonization 
and  a  mechanism  for  producing  azimuthal  asymmetries  in  the  interplanetary  dust  complex  are 
discussed:  testing  tnc  latter  may  be  possible  by  using  the  zodiacal  light  photopolarimeter  aboard 
the  International  Solar  Polar  .Mission  spacecratt  to  be  launched  in  198.1 
Suhicct  headtna.  interplanetary  medium  -  solar  system:  general  stars:  stellar  dynamics 

1.  LNJRODLCTION 

In  an  earlier  paper  tSchuerman  19721  the  effect  of  radiation  pressure  on  the  effects  e  gravitational  potential  in  the 
restricted  three-body  problem  was  investigated.  If  one  component  of  a  binary  system  is  very  luminous  (  -  l(J4  L  ).  it 
was  shown  that  the  familiar  "contact  surface"  may  not  exist  and  the  atmosphere  of  the  brighter  component  has 
access  to  a  much  larger  region  of  space  than  that  predicted  by  the  classical  Roche  potential.  In  this  study,  the 
restricted  three-body  problem,  modified  to  include  the  force  due  to  radiation,  is  investigated  with  respect 
to  the  location  and  stability  of  the  resulting  equilibrium  points.  (/  =  1  -  5).  At  these  points  the  net  force  is 
zero.  The  locations  of  the  equilibrium  points  are  determined  by  the  ratio  f/i)  of  the  force  due  to  radiation  pressure  to 
the  force  due  to  gravitation  acting  on  a  mass  m  of  radiation  pressure  cross  section  a.: 

f  Lvavdx 

O’iTcGsm' 

L .  is  the  luminosity  in  the  frequency  range  v  to  v  +  dx  of  a  star  of  mass  At.  the  ot.’er  symbols  have  their  usual 
meanings.  The  radiation  pressure  cross  section  may  be  written  AQ.  where  A  is  the  geometric  cross  section  and  Q,  is  a 
dimensionless  "efficiency  factor"  which  is  a  function  of  particle  size,  shape,  and  optical  properties.  For  a  star  of 
mass  At  and  luminosity  L  =  f  L.dx.  /?  is  roughly  determined  by  the  ratio  of  A  m  for  particle  sizes  larger  than  or  equal 
to  the  wavelength  of  the  incident  radiation.  Schwehm  ( 1976).  for  example,  shows  that  '<  is  of  the  order  unity  for 
interplanetary  particles  whose  radii  are  between  0. 1  and  I  pm. 

Spacecraft  and  artificial  satellites  can  now  be  manufactured  with  m  t  A  such  that  /if  '  I  for  these  larger 
structures.  (An  extreme  example  is  that  of  "solar  sails"  which  could  utilize  radiation  pressure  to  power  spacecratt 
through  the  solar  system. )  The  term  /?  particle  applies  here  to  either  a  small  particle  or  an  artificial  body  otuated  ne.,r 
the  Sun  or  a  star  such  that  /f  is  nonnegligible. 

The  inclusion  of  the  force  of  radiation  pressure  in  the  restricted  three-body  problem  is  presented  here  in  two  steps. 
First,  the  major  radial  components  of  the  pressure  force  are  considered  for  the  case  of  a  /f-particlc  in  the  v  icinity  of 
two.  luminous,  massive  bodies.  The  two  possibly-siable  Lagra.igian  points.  F4  and  F,.  are  shown  not  to  be  confined 
to  their  familiar  triangular  positions.  Depending  upon  the  values  of  /J.  F4  and  F<  may  be  located  anywhere  within 
the  union  of  two  circles  which  are  centered  on  the  two  stars  and  which  have  equal  rad’i  corresponding  to  the 
separation  of  the  two  stars.  Since  /J  is  a  function  of  particle  size,  so  are  the  equilibrium  positions  of  individual 
particles.  A  stability  criterion,  which  reduces  to  that  of  the  classical  problem  in  the  case  /f  — »  0.  is  derived.  In  the 
second  step,  the  Poy nting-Robertson  effect  is  introduced  for  the  case  of  one  luminous  body.  This  first-order 
rela'ivistic  effect  does  not  significantly  change  the  positions  of  L,.  but  renders  the  F4  and  F?  points  unstable  on  time 
scales  long  compared  to  the  period  of  revolution  of  the  two  massive  bodies. 

II  SOLUTION  WITH  RADIATION  PRESSURE 

In  the  restricted  three-body  problem,  a  steady  state  solution  is  sought  for  the  motion  of  a  test  particle  in  the 
gravitational  field  of  two  larger,  orbiting  bodies  of  masses  Af,  and  A fz.  It  is  assumed  that  At,  and  At :  arc  in  circular 
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orbits  about  their  common  center  of  mass  and  that  their  angular  frequency  itli  is  given  by  t2:  =  (i (A/,  »  M  i  o' 
where  a.  the  distance  between  mass  centers,  is  used  as  the  unit  ol  length  It  is  further  assumed  that  the  test  particle 
remains  in  the  piane  in  which  M ,  and  A/;  revolve.  In  the  rotating  reierence  Irame  (see  F  ig.  ]  I  whose  origin  is  the 
center  of  muss  and  w  hose  angular  frequency  is  £1.  the  two  components  ot  the  equation  ol  motion  of  a  /t-particle  with 
coordinates  A .  )  )  are 

X  2  V  (1  -  /l,)/i(A  -  u  -  I )  (1  -  /Mil  -  nll.V  -i-  ,i) 

,  '  —  A  7  ~  j  *  (  —  ) 

If-  Q  r,-  r2 

i  --V  f  (I  -  li,)u  (I  -  /Till  -  11)1.. 


\f,  is  located  at  A  •=  I  -  u.  A/  •  is  located  at  A'  =  - p.  and  it  =  A/,  I  A/,  -  A/. )  The  distances  from  the  point  ( A'.  )  I 
to  the  centers  of  mass  A/,  and  A/;  are  r,  ■=  i(A  -r  /<  -  1 );  +•  >  :  and  r:  -  it  A'  a-  u )*  -a  \  respectively.  The 

two  terms  on  the  left-hand  sides  of  equations  (2)  and  (3)  represent  the  acceleration  and  the  Coriolis  force  experienced 
be  the  />-partide  The  centrifugal  acceleration  has  been  placed  on  the  right-hand  sides  (RHSs)  ol  the  equations 
because,  like  the  remaining  forces,  it  is  a  function  of  position  only.  The  last  two  terms  on  the  RHSs  ol 
equations  (2)  and  i3i  represent  the  components  of  the  "effective"  gravitational  forces.  These  terms  difier  Irom 
those  of  the  classical  problem  in  that  each  component  has  been  reduced  hv  a  factor  ol  the  form  1  -  /)’. 

The  equilibrium  points  Li  and  L<  are  found  bv  setting  the  RHSs  ol  equations  (2 )  and  1 3 )  equal  to  zero  with  )  -  if 
The  resulting  solution  is 

r,  —  1 1  —  /q)1  '  -  0,  .  r.  =  (1  -  /M1  ’  =  <>;  .  (4 1 

or  equivalently 

I  o.  /  d,:  -  d.:  -  1  ;;T  ; 

A  -  •  -  /«  =  An .  >•=  :.Vh  1  (  ••  -  )  J  -  v  <-> 

From  solution  1 4 >  follows  that  the  /.4  and  /.<  points  exist  only  if  all  of  the  following  conditions  are  satisfied 

n,  -  A,  >  1  .  0  •;  d,  s.  1  and  0  <  d;  <  I  .  16) 

These  results  are  plotted  in  figure  2  w  Inch  shows  th  positions  of  the  I.A  and  L<  points  as  a  I  unction  ot  /f,  and  />; 
land,  equivalent K .  o,  and  <i.t.  for  <>,  -  o,  -  I .  the  classical  triangular  points  are  recovered.  For  all  other  values  ol 
A.  and  (>-  satis! v ine  conditions  16).  /.4  and  /.<  are  confined  to  the  intersections  of  the  two  sets  ol  circles  given  by 
equation  (4).  for  the  limiting  case  d.  -  — *  I.  /.4.  £.*.  and  the  interior  Lagrangian  point  L,  all  converge  to  the 

position  A'  -  1  -  o,  -  u.  )  -  (l.  For  d,  -  d,  •.  1.  /.4  and  /.<  have  no  solutions. 

The  stahiiit :  of  the  equilibrium  points  is  determined  by  investigating  the  motion  resulting  from  a  small 
displacement  ol  the  /f  part  tele  from  its  equilibrium  position.  II  equations  (2)  and  l3)  are  lineari/.ed  with  respect  to 
small  displacements  of  the  form  .V  -  A ,,( I  -  viand  )  =  >',>(  1  —  t  ).  the  resulting  linear  equations  admit  solutions 
with  an  explicit  time  dependence  of  the  form  a  =■  v  exp  (/T2/I.  v  =  y  exp  ( pdt).  Equations  (2)  and  (?)  then  become 
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Fu,  The  powiblc  positions  ot  thc  equilibrium  points  C4  fin  upper  half-plane i  and  /.«  nn  lower  half-plane  i  js  u  Inaction  ol  // ,  and  P:  or  <'  ■ 
ana  [,'  =  1 1  --  p>‘  Jj  The  heavy  line  applies  to  sys'tim  with  oni>  .me  iuminou*  body.  M. 

in  order  that  the  amplitudes  A',,.v  and  )\,  v  he  arbitrarily  small,  the  determinant  oT t his  homogeneous  set  ofequations 
must  v  anish  tie.,  equations  t  ■ )  and  (8)  must  dilfer  by  only  a  constant  factor.  This  restricts  the  possible  values  o  ip  to 
the  roots  of  the  equation 

/t1  -  p:  +  9q(]  -  ft)h  =  0  .  where  h  =  1  -  j  ---  - -  j  .  (9) 

These  four  in  =  1  -4)  roots. 

Pn  =  ±i\\  ±  [i  -  9/4(1  -  n)h]'  ;  .  (10) 

determine  the  normal  modes  of  the  linearized  system.  In  order  that  the  motion  be  stable,  none  of  the/>„‘s  can  contain 
a  positive,  teal  part  which  would  thus  lead  to  an  exp  t  /?fiM  increase  in  the  motion  of  the  /^-particle.  From  equation 
( 1 0 1  it  follows  that  the  motion  is  stable  provided  that 

1  >  36q(  1  -  q  )b  . 

Lagrange's  classical  criterion  for  stability  is  recovered  in  the  li mil  of  no  radiation  pressure  (b,  — ♦  1 .  b;  - 
for  which  equation  (11)  becomes 

1  >  27q(  1  —  q)  or  [(1  —  q)  or  q]  >  0.96148  . 

The  greatest  restriction  on  q  occurs  for  the  case  (b,2  +  b;:)— *  1  for  which  equation  (11)  reduces  to 

1  >  36q(  l-/i)  or  [( 1  —  q )  or  q]  >  0.97140  .  (13) 

Thus,  although  radiation  pressure  producesa  large  effect  on  the  location  of  the  T4  and  L?  points,  it  hasa  very  small 
effect  on  restricting  the  values  of  q  necessary  for  stability  (as  formulated  in  this  section). 

A  similar  analysis  can  be  carried  out  for  the  equilibrium  points  L,.  L2.  and  /.,  foi  which  Y  =  0  in  equations  '2) 
and  ( 3 ).  Itis  found  that,  as  in  the  classical  case,  these  points  arc  unstable  and  no  equilibrium  points  exist  for  >  I 
or  /I.  >  1  The  sole  exception  to  this  result  is  the  degenerate  case  for  which  (b,  ■+•  b:)  — *  1  :  then  L,  the  interior 
Lagrangian  point  and  L4  and  L5  all  converge  to  the  same  point  on  the  A  axis  ( I  -  b,  -  n ).  and  that  point  is  stable 
lor  all  values  of  q. 

Itis  of  interest  to  consider  the  implications  of  these  results  for  the  solar  system.  The  subscript  1  then  refers  to  the 
Sun  while  subscript  2  refers  to  a  planet.  In  this  case  j!:  vanishes,  b ,  is  unity,  and  the  subscripts  on  /},  and  «5,  are 
dropped  for  convenience  (here  and  throughout  the  rest  of  this  paper).  Submicron  particles  have  values  of  [i  which 
can  he  an  appreciable  fraction  of  unity.  Such  particles  c  annot  be  "trapped"  at  the  classical  triangular  points  as  are 
the  Trojan  planets  of  the  Sun-Jupiter  system.  Even  if  other  forces  (e.g..  Lorentz.  Poyr.ting-Robertson)  are  not 
considered,  the  classical  triangular  points  are  not  eouilibrium  points  for  these  small  particles.  Their  equilibrium 
points  he  on  the  arc  of  a  circle  w  hich  is  centered  on  A/,  and  w  hose  radius  is  the  d.  stance  separating  the  Sun  and  M:. 
This  arc  is  show  n  as  the  heavy  curve  in  Figure  2.  The  location  of  the  equilibrium  points  on  this  arc  are  parametrized 
by  /i  As  /(  — >  1  the  /.4  and  /.<  point*  approach  the  center  of  the  Sun  If  no  other  forces  are  considered,  this 
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like  .1  si/c  or  mass  spectrometer .  a  distribution  ol  -ni.ii!  partuie -ura: 
cure  .1  would  be  maintained  tor  an  arbitrarily,  lone  '.me  lit  '.tie  ;re\ 
to  one  body,  is  taken  into  account  When  tins  toiee  is  coiisidem:  :! 
i'll  time  scales  lone  compared  to  il  1 

aned  thus  tar  have  an  interestin'.:  application  lot  artificial  '.i  c 
■en  suggested  (e  e  .  O'Neill  I V74 1  that  the  classical  tr  i.ineuiai  points  oi 
be  convenient  sites  to  locate  future  space  colonies  Such  coiomes 
least  some  of  their  raw  materials  Present  das  technology,  permits  the  c. 

!'  can  be  made  or  varied  in  the  range  I)  <  ii  <  1  i.A  spacecraft  with  /;  ■  1 
::mg  "  vehicle  lor  chasing  Halley. N  comet  [  I  riedman  ll>~6]  i  A  space  co 
oui.i  thus  "park  1 1  e. .  t  emain  si  at  lonarv  with  respect  to  t  he  Sun  amt  a  p 
:  energy.  >  not  only,  at  the  classical  positions  but  at  an i  heliocentric  distance 
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P. n.'.mg  i  ! -a < '  w as  the  nrst  to  suggest  that  /{-particles  in  solai  orbit  sutler  a  gradual  loss  ol  angular  momentum 
,i::C  ultimate!',  spiral  into  the  Sun  In  a  svstem  ol  coordinates  where  the  Sun  is  at  rest,  radiation  scattered  by.  the  /)'- 
particle  in  me  sines".: on  .'!  motion  sutlers  a  blueshift.  and  in  the  opposite  direction  it  is  redsht! ted  I  his  gives  rise  to  a 
r.e:  drag  force  which  opposes  the  direction  ol  motion.  The  proper  relativistic  treatment  ol  this  problem  was 
!  o  i  in  a  .a  i  cd  by.  R.-bei'.son  d'T'D  who  showed  that  to  tirst  order  in  I  c.  the  radiation  pressure  torce  is  given  by. 


/.(/A/  [  r.  I  f  •  r.  lr.  I 
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(14) 


I  i.s'  .vr.  Tdc:  has  been  treated  :n  the  preceding  sections.  The  second  term  on  the  right-hand  side  of  equation 

!4i  represents  due  Doppler  shut  ol  the  incident  radiation,  and  the  last  term  is  the  Povnnng  drag.  These  last  two 
:r. -  'aken  together  are  tiie  choc:  under  consideration. 

I  ••  mcorpota'.e  the  Powuing-Rohertson  force  into  the  restricted  three-bods  problem  as  formulated  in  the 
preceding  'wivi:  the  semens  vector  l  'of  the  />'  particle  must  be  traristormed  to  the  rotating  s> stem  centered  on  the 


I  -  v  -  ii  ■  R  .  v  —  A i  -  Yj  .  (15) 

w  here  i  and  j  me  ,.::it  or; hogonal  v eetors  m  the  rotating  frame  ( see  Tig.  1 ).  In  dimensionless  units  the  /./components 
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where  .  which  .ontams  the  parameter  analogous  to  ID.  is  given  by.  y  s.  jUiaL 2  c  The  positions  of  the  La  and  /.. 
pi  i::;.  arc  once  ucain.  found  by.  setting  the  RHSs  of  equations  ilbi  and  11')  equal  to  zero  with  )'  *  0  To 

Mr  s’  ,  rue:  it:  .  1  he  s.  lulloii  is 
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\  an.!  )  rep:  e  cm  I  he  eqmli  Irr  mm  positions  of  the  and  /.^  points  w  hen  the  I’ov  mine- Robert  son  e  licet  is  noi 
in!  i  led  and  w  hen  /;_  ♦  0.  n,  -I  in  equations  (5): 

A  ,,  I  //  -  <>-  2  .  T(1  y  M  1  4 11  r  *  (1  (21  ) 

1  sjsiiiun;,  '  Is  i  id'll  are  ..did  provided  the  first  order  terms  in  ;  are  small  compared  to  units 

I  siai d;tv  ■  A  th-  .leads  stale  solution  derived  above  is  analyzed  bv  the  same  technique  as  that  used  in  A  11 
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Equuuons  <  16)  and  (17)  are  linearized  with  respect  t ii  small  displacements  about  the  equilibrium  points.  and  . 
given  by  equation  (19)  and  (20): 

A' =  A0  [I  -r  ,v  exp  t /nzrl]  .  }' =  }„[l  -r  t  exp  (pO/)]  .  (22) 

A  straightlorward  calculation  yields  the  result,  analogous  to  equations  (7)  and  (8).  that 


■Vo". -v‘,r 
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For  a  solution  to  exist,  the  determinant  of  this  set  of  linear,  homogeneous  equations  must  vanish.  This  yields  a 
fourth  order  equation  in  p  of  the  form 

(1  +  -f  a  ,/r  -HI  -  a:)p:  -f-  a,p  +  9p<  1  -  p)b(  1  -  u,,|  =  0  .  (25) 

w  here  h.  defined  in  equation  (9).  reduces  to  b  =  1  -  <r  4  and  at  (j  —  0.  1 . 4 1  must  be  first  order  terms  in  y.  We 

seek  a  solution  of  equation  (25)  of  the  form  p  =  p„(  1  -f  x  +  ip)  where  p„  is  given  by  equations  ( 10).  restrictions  (II) 
applies,  and  both  x  and  p  must  be  first  order  terms  in  y.  Then  p  takes  the  form 

p  =  ~i(\  -r  2)[v  ±({  ~  9.11(1  -  p)h)'  ']'  2  +  i;[T  +  (i  -  Q;i<!  -  p) b)‘  2  (26) 

Thus,  if  the  /;’  particle  is  displaced  from  its  equilibrium  position,  it  will  oscillate  w  ith  frequencies  determined  by  the 
imaginary  part  of p  in  equation  (26)  while  the  grow  th  m  the  amplitude  of  those  oscillations  is  determined  by  the  real 
part  of  p.  It  is  the  time  scale  of  the  latter  that  will  be  investigated  here.  If  equation  (26)  is  substituted  into  equation 
(25).  the  following  expression  for  p  is  obtained: 

[*  ±  (1  ~  9p(l  -  g)b)'  •]■  ±  (j  -  9p(  1  -  p)h)'  -  a,l  p?) 

n  T  ±4(i  -  9/i(l  -  p)b)1  2[f  ±  (i  -  9/t(l  -  p)b)'  2] 

To  evaluate  p.  it  is  only  required  that  a}  and  a,  ( the  coefficients  of  p3  and  p)  be  computed  rather  than  the  entire 
determinant  of  the  set  of  equations  (23)  and  (24).  To  first  order  in  y.  wc  find 

3y  3y  f ,  ,  >o2(l  ~p)"| 

Cl}  =  —  •  0,  =  -  —  1-1 - 77 -  •  (28) 

(f  d-  |_  62 

so  that  the  e-folding  time  (  T)  for  the  growth  of  /1-particle  oscillations  about  the  equilibrium  points  T4  and  Ls  is  given 
bv 


1  yfi  p-jj  ±  [»  -  9.n(l  -  ,<)(!  -  S:  4)]1  2j  +  3[1  -Hi  -  /()( 1  - 

T~  C  (P)  ~  4b 2  ±[i  -  9p(l  -  p)(l  -  b2/4)]'  : 


ZJ)T 


Under  restriction  (11).  T  is  always  positive  for  at  least  one  choice  of  mode  (sign)  in  equation  (29).  Therefore,  / i - 
particles  in  general  have  no  stable  equilibrium  solutions  of  the  restricted  three-body  problem  if  the  Poynting- 
Robertson  elfect  is  taken  into  consideration. 

T  takes  a  more  convenient  form  when  applied  to  the  solar  system.  For  all  of  the  planets  including  Jupiter,  terms 
which  contain  the  factor  (1  —  p)  can  be  ignored  in  equation  (29).  The  mode  with  the  smallest  characteristic  time  (the 
fastest  growing  oscillation)  has 


[\~P)2\a2  (\  -  P)2 


544u2  yr  (a  in  AU)  . 


3Qy  3pMeG  P 

A  convenient  and  representative  value  for  p  of  0.57  makes  (1  -  P)2  i  'P  =  1.  Thus,  for  example,  the  Sun-Jupiter 
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v\ stent  t a  =  5.2  -\U  I  results  in  7  i  14. "IMI  yr  This  meats  that  ti  n  particle  with  //  =  <152  isdisplaced  slightly  from  its 
equilibrium  position,  it  will  oscillate  about  that  position  in  some  linear  combination  ol  its  tour  normal  modes  and 
iIk  amplitude  oi  that  motion  will  increase  by  a  factor  off  in  I4.~0()\r 


tv,  disi  t  sstos 

It  ;s  not  at  all  evident  that  the  above  mechanism  could  lead  to  a  significant  increase  in  the  density  of  interplanetary 
//-particles  along  the  arc  schematically  shown  as  a  heavy  line  in  Figure  2.  Other  forces  act  on  these  //-particles  which 
cannot  be  treated  by  the  analytical  methods  employed  here.  Interactions  with  the  gravitational  fields  of  intervening 
planets  could  displace  some  of  the  //-particles  far  enough  from  their  equilibrium  positions  that  the  Coriolis  force  : 
insufficient  to  restore  them.  This  effect  would  leave  large  gaps  in  the  arc  ot  //-particles.  Jupiter,  because  of  its  great 
mass,  would  be  particularly  effective  in  disrupting  the  dust  in  arcs  associated  with  the  other  planets ;  conversely  .dust 
in  the  Sun-Jupiter  arc  would  be  the  least  affected  by  planetary  encounters.  The  Lorentz  force  due  to  the  solar 
magnetic  held  also  acts  on  (charged)  //-particles.  This  conservative  force  is  roughly  40'^,  that  of  the  solar  gravitation 
(Greenberg  and  Schuerrrun  1975/ ) .  The  uncertainty  in  the  charges  of  the  //-particles  and  in  the  long  term 
fluctuations  and  reversals  of  the  solar  magnetic  field  at  distances  greater  than  1  AU  make  the  inclusion  ot  the 
Lorentz  force  in  this  problem  particularly  difficult.  Finally,  the  ions  in  the  solar  wind  impinging  on  the  //-particles 
and  the  Coulomb  interaction  between  solar  wind  electrons  and  the  //-particles  contribute  nonconservative  drag 
forces  which  together  are  estimated  to  be  about  30"”o  of  the  Poynting- Robertson  effect  ( Bandcrmann  1 467 .  Mi  scorn 
1976«).  These  forces  must  surely  operate  to  shorten  the  time  scales  (T)  in  equation  (30). 

Even  if  the  planetary  encounters  do  not  disrupt  the  dust  along  the  entire  arc.  and  even  if  the  large  but  conservative 
Lorentz  force  averages  out  to  a  small  effect  over  manv  solar  cycles,  there  still  remains  the  need  for  a  replenishment 
mechanism  if  an  increase  in  the  density  of  //-particles  along  a  planetary  arc  is  to  be  expected  on  a  "steady  state" 
basis.  The  major  source  of  replenishment  for  the  interplanetary  dust  complex  is  thought  to  be  cometary  debris 
(Whipple  1955.  1967;  Dohnanyi  1970.  1972.  1973),  Comets  shed  both  //-particles  (some  with  //  >  1 )  and  larger 
particles;  light  scattering  by  the  latter  is  the  main  contribution  to  the  zodiacal  light.  These  larger  particles  ultimately 
fragment  through  self-collisions  (Zook  1975)  or  possibly  by  rotational  bursting  (Paddack  and  Rhee  1975;  Misconi 
19766).  It  is  not  evident  whether  or  not  these  sources  of  //-particles  provide  a  wide  enough  distribution  in  real  space 
and  velocity  space  to  deposit  a  significant  fraction  of  "zero- velocity  "  //-particles  at  their  appropriate  equilibrium 
positions  in  times  short  compared  to  T  in  equation  (30). 

Mechanisms  which  tend  to  produce  azimuthal  asymmetries  in  the  interplanetary  dust  complex  are  of  great 
interest  because  present  methods  of  inverting  spacecraft  observations  of  the  zodiacal  light  (Dumont  1973; 
Schuerman  1979)  depend  upon  the  assumption  of  axial  symmetry  for  the  dust  cloud.  Fortunately,  the  International 
Solar  Polar  Mission  whose  launch  is  scheduled  for  1983  February  provides  the  opportunity  to  assess  this  assumption 
directly  .  A  photopolarimeter  (jointly  designed  by  teams  at  the  Ruhr-University.  Bochum.  FRG.  and  the  State 
University  of  New  York  at  Albany )  will  map  the  brightness  and  polarization  of  the  zodiacal  light  throughout  the  4) 
year  mission  to  Jupiter  and  back  over  (under)  the  poles  of  the  Sun.  During  the  two  periods  of  solar  polar  passage,  the 
spacecraft  will  be  in  the  center  of  the  solar  system  but  at  1.7  AU  above  (and  below)  the  ecliptic  plane.  From  these 
vantage  points,  isophotes  of  the  zodiacal  light  will  appear  as  circles  only  if  the  interplanetary  dust  complex  is  truly 
azimuthaliy  symmetric.  In  particular,  gross  enhancements  in  the  spatial  density  of  the  dust  along  arcs  associated 
with  planets  should  be  detectable.  Furthermore,  this  study  predicts  that  any  diffuse  light  associated  with  the 
classical  libration  (Z.4  and  Z.5)  points  ol  the  pianei-iun  systems  must  be  due  to  grains  larger  than  about  a  micron 
because  the  triangular  points  are  not  equilibrium  points  for  //-particles. 

This  w  ork  was  supported  by  the  Air  Force  Office  of  Scientific  Research  (contract  F49620-78-C -0013)  and  by  the 
National  Aeronautics  and  Space  Administration  tJPL  contract  955137). 


RELERLNCES 


Handermann.  L  W  ]967.  Ph  D  thesis.  1.  niversity  ol  Maryland 
Dohnanyi.  J  S.  1970.  J.  Geophys.  Res..  75.  3648 
1972.  Icarus.  17.  1. 

------  1973.  in  Evolutionary  and  Physical  Properties  ot  Meteoroids. 

ed  Hemenwau  Miliman.  and  Cooks  (Washington :  NASA  SP 
3 1 9  >.  p  .363. 

Dumont.  R  1973.  Planet  Spate  Set..  21.  2149 
Friedman,  L  D  1975.  in  Proceedings  of  the  Shuttle  Based  ( Dmetarx 
Science  Workshop,  ed  (j  A  Gary  and  K  S  Chiton  (Huntsville 
NASA  Marshal/  Space  Plight  Center  1.  p.  251. 

Greenberg.  J  M  .  and  Schuerman.  D  W  1978.  Suture.  275.  39. 
Misconi.  N  V  1976ij.  Astr.  Ap..  51.  357 

19766.  (/eophvs  Res.  letters.  3.  585. 


O'Neill.  G  k.  1974.  Plu  s  Todas .  27.  32 

Paddack.  S  J..  and  Rhee.  J  W  19  75.  Geophys.  Res  Letters.  2.  565 
Poynting.  J  H  1905.  Phil  Trans  Row  Soe  .  A.  202.  525. 
Robertson.  H  P  1937.  V/  SR  A S.  97.  423 
Schuerman.  D  W  1972.  Ap  Space  Set  .  19.  351 
- 1979.  Planet.  Space  Set..  27.  551. 

Schwehm.  G  1975.  in  Lecture  Soles  in  Phvstcs.  Vol  48.  ed  H 
Eisasser  and  H.  Fechtic  (Heidelberg  Springer- Verlag).  p  459 
Whipple,  F  L  195 5.  Ap  J..  121  75Cf 

— -  1967.  in  The  Zodiacal  Liahl  and  the  Inter  planet  arc  Medium. 

ed.  J  L.  Weinberg  (Washington  NASA  SP  150).  p  409 
Zook.  H  A  1975.  Planet  Space  Set  .  23.  1391 


Dos all>  W.  Sem  lrmas:  Space  Astronomv  Laboratorv.  Slate  University  of  New  York  at  Albanv,  Executive  Park 
East.  Albany.  NY  12203 


THE  PHOTOMETRIC  CENTER  OF  THE  GEGENSCHEIN 


1  . 


by 


Nebil  Y.  Misconi 


Space  Astronomy  Laboratory 
University  of  Florida 
1810  N.W.  Sixth  Street 
Gainesville,  Florida  32601 


Seven  Manuscript  Pages 
Five  Figures 
No  Tables 


3. 


j  abstract 

Model  calculations  are  used  to  evaluate  two  factors  which  de- 

ternine  the  position  of  the  photometric  center  of  the  Gegenscheinr 

the  increased  scattering  efficiency  of  the  interplanetary  dust  near 

o  c 

backscattering  (scattering  angles  6  'v  165  -  180  ),  and  the  spatial 

density  distribution  of  the  dust.  Computer-generated  brightness 
contours  are  used  to  investigate  which  of  these  two  factors  domi¬ 
nates,  using  empirical  scattering  functions  with  and  without  a 
brightness  enhancement  in  the  backscattering  region.  It  is  found 
that  the  effect  of  the  enhanced  scattering  of  light  by  dust  in  the 
backscattering  region  overrides  the  effect  of  the  spatial-density 
distribution  of  the  dust.  As  a  result,  the  photometric  center 
should  be  observed  at  the  antisolar  point  at  all  times,  except 
possibly  when  the  antisolar  point  is  at  its  maximum  displacement 
from  the  symmetry  plane. 


1.  INTRODUCTION 


k  . 


The  position  of  the  photometric  center  (maximum  brightness) 
of  the  Gegenschein  has  been  thought  to  relate  directly  to  the 
spatial  density  distribution  of  interplanetary  dust  in  the  re¬ 
gion  of  the  Gegenschein.  From  observations  of  the  Gegenschein, 
Dumont  (1965)',  Robley  (1965),  and  Roosen  (1968,  1970a,  1970b) 
found  the  photometric  center  to  be  close  to  the  ecliptic  plane. 
However,  Tanabe  (1964,  1965,  1980)  and  Wclstencroft  and  Rose  (1967) 
found  the  photometric  center  displaced  from  the  ecliptic  plane 
and  close  to  the  invariable  plane  of  the  solar  system.  Earlier 
observations  (1950's)  shewed  similar  displacements  (Figure  1). 
Pioneer  10  and  11  observations  (Weinberg  et  al . ,  1973)  showed  be¬ 
yond  doubt  that  the  Gegenschein  arises  from  enhancement  in  the 
backscatter ing  of  interplanetary  dust.  The  question  then  arises 
as  to  whether  the  position  of  the  photometric  center  is  directly 
related  to  the  maximum  dust  density  distribution.  We  would  like 
to  emphasize  that  there  are  two  factors  which  determine  the  posi¬ 
tion  of  the  photometric  center:  the  increased  scattering  effi¬ 
ciency  of  the  interplanetary  dust  near  backscattering  (scattering 
angles  of  approximately  165°  to  180°)  and  the  spatial  density  dis¬ 
tribution  of  the  dust.  In  other  words,  can  the  position  of  maxi¬ 
mum  brightness  be  at  the  antisolar  point  if  the  plane  of  maximum 
dust  density  does  not  contain  this  point?  Computer-generated 
brightness  contours  are  used  to  investigate  this  question  using 
empirical  scattering  functions  with  and  without  a  brightness  en- 


hancement  in  the  backscatter ing  region,  thereby  separating  the 
effect  of  the  increased  brightness  from  the  effect  of  the  spatial 
density  distribution.  It  will  be  shown  later  that  the  enhance¬ 
ment  in  the  hackscattering  overrides  the  spatial  density  distri¬ 
bution  of  the  dust  in  the  region  of  the  Gegenschein. 


0 


MODEL  CALCULATION 


Predictions  were  made  earlier  (Figure  1)  of  the  displacement 
in  the  ecliptic  latitude  of  the  photometric  center  for  the  plane 
of  maximum  dust  density  (symmetry  plane)  at  the  invariable  plane 
or  at  the  orbital  plane  of  Mars  (Misconi,  1977).  These  predic¬ 
tions  were  based  only  on  the  effect  of  the  spatial  density  distri¬ 
bution  of  the  dust.  Using  ground-based  observations  of  the  Gegen- 
schein  (Figure  2)  Dumont  (1965)  found  the  photometric  center  of 
the  Gegenschein  to  be  in  the  ecliptic  plane  and  the  photometric 
axis  (locus  of  points  of  maximum  brightness)  to  be  near  the  in¬ 
variable  plane.  Dumont's  findings  agree  with  our  predictions  at 
15°  east  and  west  of  the  antisolar  point  (Misconi,  1977),  but 
disagree  with  the  predictions  in  Figure  1  (where  the  effect  of  en¬ 
hanced  backscattering  of  the  dust  was  not  evaluated).  Using  a 
computer  code  developed  earlier  by  Hanner  (Misconi  and  Hanner, 

1975),  an  attempt  was  made  to  reproduce  Dumont's  observational  map 
of  the  Gegenschein  (Figure  3).  This  map  was  generated  using  Dumont* 
(1975)  empirical  scattering  function,  which  he  deduced  from  his 
Gegenschein  brightness  data. 

To  separate  the  effect  of  increased  scattering  efficiency  of 
interplanetary  dust  near  backscattering  from  the  spatial  density 
distribution  of  the  dust,  a  computer-generated  brightness  contour 
map  was  made  (Figure  A)  using  the  empirical  scattering  function 
without  a  brightness  enhancement  in  the  backscattering  region. 

The  spatial  density  distribution  above  and  below  the  ecliptic 


7  . 


plane  was  assumed  to  follow  the  fan  model  suggested  by  b^.i..ert 
ei  cl.  (1976).  The  fan  model  was  intended  to  serve  as  a  good 
example,  although  there  are  several  other  models  that  have  also 
been  suggested  (see,  for  example,  Geise,  1972,  and  Leinert  et  al . , 
1976).  Choosing  a  specific  model  for  the  dt  .sity  distribution  is 
not  considered  critical  for  the  purpose  of  separating  the  two  ef¬ 
fects  mentioned  above.  This  follows  from  Figure  1,  since  the 
displacement  in  ecliptic  latitude  of  the  symmetry  plane  is  small 
(2.5°)  for  the  time  of  Dumont's  observations  (February  9),  and 
the  various  models  of  the  dust  distribution  do  not  differ  signi¬ 
ficantly  at  small  latitudes  from  the  ecliptic.  Figure  4  shows 
very  clearly  that  if  it  were  not  for  the  effect  of  enhanced  back- 
scatter,  the  photometric  center  of  the  Gegenschein  would  be  dis¬ 
placed  from  the  antisolar  point  by  nearly  the  same  displacement 
of  the  symmctiy  plane. 

The  question  now  arises  as  to  whether  or  not  the  photometric 
center  continues  to  be  seen  at  the  antisolar  point  at  all  seasons, 
even  during  the  maximum  displacement  of  the  symmetry  plane  from 
the  ecliptic.  As  the  displacement  of  the  symmetry  plane  from  the 
ecliptic  increases,  the  decrease  in  the  number  density  of  the  dust 
at  the  antisolar  point  becomes  model  dependent.  Therefore,  the 
answer  in  this  case  depends  heavily  on  the  actual  spatial  density 
distribution  of  the  dust  above  and  below  the  ecliptic  plane.  A 
computer-generated  map  (Figure  5)  was  made  using  the  fan  model  men¬ 


tioned  earlier  for  a  time  when  the  antisolar  point  is  at  maximum 
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displacement  in  ecliptic  latitude  from  the  invariable  plane  (this 
corresponds  to  A^  'v  350°  -  50°  and  A  'v  170°  -  220°)  i.e.,  during 
March,  April,  September,  and  October.  Figure  5  shows  that  if 
one  uses  the  fan  model  the  photometric  center  will  not  be  seen  at 
the  antisolar  point;  it  will  tend  instead  toward  the  same  displace¬ 
ment  from  the  ecliptic  as  the  symmetry  plane  near  the  antisolar 
point.  This  also  means  that  the  effect  of  enhanced  backscatter ing 


is  no  longer  apparent. 


3.  DISCUSSION  AND  CONCLUSION 


This  study  addresses  the  present  confusion  in  published  re¬ 
sults  on  the  position  of  the  photometric  center  of  the  Gegenschein. 

The  significance  of  the  true  position  to  our  understanding  of  the 
dynamics  of  interplanetary  dust  has  been  discussed  earlier  (Misconi, 

1977).  It  is  concluded  from  our  model  calculations  that:  (a)  The 
effect  of  the  enhanced  scattering  of  light  by  dust  in  the  back- 
scattering  region  (6  ^  165°  -  180°)  overrides  the  effect  of  the 
spatial-density  distribution  of  the  dust  in  this  region;  (b)  Al¬ 
though  the  photometric  center  is  observed  near  the  antisolar  point, 
this  does  not  mean  that  the  dust  density  reaches  a  maximum  at  that 
point,  but  rather  that  the  density  has  a  maximum  at  the  symmetry 
plane;  (c)  As  a  result  of  (a) ,  the  photometric  center  should  be 
observed  at  the  antisolar  point  at  all  times,  except  possibly  when 
the  antisolar  point  is  at  maximum  displacement  from  the  symmetry 
plane.  This  possibility  depends  heavily  on  the  dust  distribu¬ 
tion  above  and  below  the  symmetry  plane.  Information  on  the  ac¬ 
tual  nature  of  the  dust  distribution  will  be  obtained  by  the 
zodiacal  light/background  starlight  experiment  using  observations 
obtained  during  the  International  Solar  Polar  Mission  (Schwehm  et  al .  ,  1981) 
Such  information  can  verify  or  deny  the  possibility  mentioned  in 
(c) .  We  also  hope  to  obtain  new  brightness  mans  of  the  Gegen¬ 
schein  at  different  times  of  the  year  from  Mt.  Haleakala,  Hawaii, 
to  establish  the  relative  positions  of  the  photometric  center  and 


the  symmetry  plane  of  the  Gegenschein  at  all  seasons. 


10. 


If  future  ground-based  observations  of  the  Gegenschein  show 
the  symmetry  plane  to  be*  near  the  invariable  plane  at  all  seasons 
with  the  photometric  center  always  at  the  antisolar  point  (in  the 
ecliptic  plane) ,  then  it  is  very  likely  that  the  dust  seen  in  the 
region  of  the  Gegenschein  is  affected  gravitationally  by  Jupiter. 


ACKNOWLEDGEMENTS 


I  would  like  to  thank  Dr.  J.  L.  Weinberg  for  critically  review¬ 
ing  the  manuscript.  This  research  was  supported  by  the  Air  Force  Office 
of  Scientific  Research  end  the  National  Aeronautics  and  Space  Administra¬ 
tion  (Planetary  Atmospheres  Program). 


Dumont,  R.  (  1965).  Ann.  Astrophys .  2£_,  263-320. 

_ .  (197  5)  .  Astror..  Astropr.ys.  IS,  405-412. 

Elsasser,  H.  ,  and  Siedcntopf ,  H.  (  1957).  Z.  Astrovhys .  41 ,  2,  132-143. 

Giese,  R.  li.  (  1972).  Zodiacal  light  and  interplanetary  particle  num¬ 
ber  densities.  In  Space  Research  XII ,  up.  437-443.  Akademie-Ver lag , 
Berlin. 

Gindilis,  L.  H.  (I960).  Soviet  Actron.  AJ .  £,  1004-1006. 

Hoffr.eister ,  C.  (  1940).  As  iron.  Sac  nr.  271  ,  49-67. 

Leinert,  C.,  Link,  H.,  and  Pitz,  E.  (  1976).  Astror..  Astrophys.  47. 
221-230. 

Lillie,  C.  F.  (1968).  An  empirical  determination  of  the  interstellar 
radiation  field.  Ph.D.  Thesis,  University  of  Wisconsin. 

Misconi,  N’.  Y.  (  1977).  Astror..  Astrophys.  61,  497-504. 

_ ,  and  Hanner,  M.  F.  (  1975).  Planetary  Space  Sci.  23 , 

1329-1335. 

Roach,  R.  E.,  and  Rees,  M.  H.  (1956).  The  absolute  photometry  of  the 
Gegenschein .  In  The  AirgtoV  and  the  Aurorae  (E.  B.  Armstrong  and 
A.  Daigarno,  Eds.),  pp .  142-155.  Pergaraon  Press,  London  and  New  York. 

Rohley,  R.  (  1965).  Ann.  Geophus.  22_,  505-513. 

Reuser. ,  R.  G.  (1968)  .  Icarus.  9,  429. 

_ .  (  1970a).  The  Gegenschein.  Ph.D.  Thesis,  University  of 

Texas. 


.  (1970b).  Icarus.  1 


184-201 . 


1 


3 . 


c  . 


Schwehr.,  0  .  H.  ,  Giese,  R.  H.  ,  Giovane,  F.  ,  Schuerir.an,  D.  V.'.,  and 
Vein'ror.c,  .1.  L.  (In  press).  Recent  developments  in  space-borne 
zodiacal  1  ight  photometry.  In  Zvs.SC  Research . 

ai.abe,  H .  ( i  .  Photoelectric  observations  of  the  Geitenschein. 
Ph.l:.  diesis,  University  of  Tokyo . 

.  (1965).  Ptti  * .  Asiron.  Zee.  Zaven.  e  ,  339-366. 

Takec'ni,  A.,  and  Miyashita,  A.  (1980).  Photometric  axis 
measurements  of  the  zodiacal  light  at  large  elongations.  In  IAV 

.Up.  PC,  Sc Z id  Particles  in  the  Solar  Easier'.  (I.  Halliday 
5.  A.  McIntosh,  Eds.),  pp .  45-48 .  Reidel  Publishing  Company, 

!  r  .  h  t  ,  Hoi  land  . 

'/•inhere,  .  ■ .  ]. .  ,  banner,  M.  S.,  Mann,  H.  M.,  Hutchison,  P.  B.  ,  and 
Fimmcl,  R .  (1973).  Observations  of  zodiacal  light  from  the  Pioneer 
10  Asteroid-Jupiter  probe.  In  Space  research  XIII  (M.  J.  Rycroft 
and  S.  K.  Runcorn,  Eds.),  pp.  1187-1193.  Akademie-Ver lag ,  Berlin. 
Woistcncrof t ,  R.  P  ,  and  Rose.  L.  J.  (  1967).  As  trophy c.  J.,  147 , 


271-292. 


FIGURE  CAPTIONS 


Fig.  1.  Predicted  positions  of  the  photometric  center  of  the 
Gegenschein,  for  planes  of  symmetry  corresponding  to  the  orbital 

plane  of  Mars  ( - )  and  the  invariable  plane  ( - )  .  Observed 

positions  of  the  photometric  center  are  also  shown. 

Fig.  2.  Isophote  map  of  the  Gegenschein  obtained  by  Dumont  at 
Tenerife  on  February  9,  1969  (Dumont,  1965). 

Fig.  3.  Computer-generated  isophote  map  of  the  Gegenschein  us¬ 
ing  Dumont's  empirical  scattering  function  (Dumont,  1975). 

Fig.  9.  Brightness  contours  of  the  Gegenschein  generated  with 
no  enhancement  in  the  backscattering  region  (0  n,  165°  -  180°). 

Fig.  5.  Computer-generated  brightness  contours  for  times  when 
the  symmetry  plane  is  near  its  maximum  displacement  from  the  eclip¬ 
tic  plane. 
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SHORT  COMMUNICATION 


ON  THE  EQUIVALENCE  OF  EXPRESSIONS  FOR  INVERTING  THE 
3*  G*in£.NSIONAL  ZODIACAL  LIGHT  BRIGHTNESS  INTEGRAL 


Dona-ld  W.  Schue man 


Space  Astronomy  Laboratory 
State  University  of  Nev  York  at  Albany 
Executive  Park  East 
Albany,  Nev  York  12203 

(deceived  9  September  1980) 


Abstract  -  In  the  research  notes  of  this  Journal,  both  Schuerman  ( 1979 , 
henceforth  Paper  I)  and  Buitrago  ( 1979  ,  henceforth  Paper  II)  independently 
derived  expressions  for  the  general^,  mathematical  inversion  of  the 
zodiacal  light  brightness  integral.  In  this  communication,  it  is  shown 
that  the  expressions  are  equivalent,  differing  only  in  the  choice  of 
reference  systems. 


In  Paper  I,  a  system  of  coordinates  (r  ,h  ,X*,8)  is  employed  where  (r  ,h  )  designates  the 
r  0*0*0*  0*0 

position  of  the  observer  in  cylindrical  ecliptic  coordinates  and  (X^,0)  defines  the  "Zook" 

direction  in  sun  centered  (X  =X-X  )  ecliptic  coordinates.  It  was  in  this  reference  system  that 

o  o 

the  3-dimensional  inversion  of  the  zodiacal  light  brightness  integral  was  first  formulated 

(Schuerman,  1977).  In  Paper  II,  Buitrago  used  a  different  set  of  independent  variables  (F,0Q,e,i) 

to  derive  a  result  which  was  difficult  to  compare  with  the  brief,  published  abstract  of  Schuerman 

(1977);  Paper  I  was  in  press.  In  Paper  II,  the  position  of  the  observer  is  designated  by  (R,Bq) 

where  R  is  the  observer’s  heliocentric  distance  and  B  is  his  ecliptic  latitude.  The  observer’s 

o 

J  "look"  direction  is  defined  by  c,  the  elongation  angle,  and  i,  the  inclination  of  the  scattering 
|  plane  (see  diagrams  in  Papers  I  and  II/. 

|  In  both  papers  it  is  sought  to  invert  the  brightness  integral  (called  Z  in  Paper  l) 


R2 

-v 


no  dl  , 


(1) 


1 


The  only  assumption  being  that  the  interplanetary  dust  complex  is  azimuthally  symmetric. 
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vhere  F  R*/F*  is  the  solar  flux  density  (ergs  cafes'*)  incident  upon  a  dust  number  density 
n  (cm~3)  whose  average  differential  scattering  cross  section  is  o  (cm2sr~l).  The  integral  is 
over  the  line  of  sight,  dl.  In  order  to  simplify  the  notation,  ve  vill  use  i>  instead  of 
Buitrago's  B^,  and  ve  drop  the  subscripts  and  the  prime  on  the  author's  independent  variables: 

* ,  3*0"*  r,  h  h,  X^**  X.  If  the  integrand  of  (l)  i6  designated  by  £,  the  brightness /unit 
line  of  sight,  the  results  of  both  papers  can  be  summarized  as  follows: 


|  Paper  11 
l  eq.  23 J 

dr,h,A,B) 

31  31 

*  - —  cosB  cos X  -  —  sine  - 
fr  3b 

31 

— rr  cosB 
rdX 

sinX  , 

(2) 

[Paper  II) 
l  eq.  5  1 

C!  ,i) 

3i  3i  ai 

c  Tr  COSE  -  — -  sine  -  ~T 

dp  Hdt  Rd^ 

sine  ^1- 

cos2l  y* 
c  os‘W 

(3) 

It  is  the  purpose  of  this  cosssuni  cation  to  point  out  that  (2)  and  (3)  are  equivalent.  The 
two  sets  of  independent  variables  are  related  by 


*  h/r  , 

(M 

ccsf  CCS  A  , 

(5) 

(6) 

Equations  (U)  are  obvious  from  the  coordinate  '  •:  *:cr.s;  (5)  and  (6)  are  derived  in  Paper  X. 
The  trans  formations  of  the  partial  derivatives  < '■  ate  simplified  by  noting  that  for  all 

elements  along  the  li ne-of-sight ,  both  0  and  i  are  constants;  thus  3I/3B  *  3I/3i  =  0.  The  chain 
rule  for  the  parti als  in  (2)  is  then  given  by 


3:  31  3?  3:  3i  t  31  3c 

3  r  or  or  3  *  3  r  3  c  3  r 

3 1  3 1  3_F  3T  3 %,  31  3c 

3h  3 p  3t  3*  3h  3c  3h 


21 

3  a 


0 


0 


3c 
3  X 


(7) 


The  null  elements  reflect  the  fact  that  3H/3X  *  3^/dX  =  0;  i.e.,  the  position  coordinates  are  in¬ 
dependent  of  the  "looX*  coordinates. 

The  transformation  of  (2)  is  nov  straightf orvard.  As  an  example,  ve  compute  the  coefficient 
of  3I/3i^  in  (3).  Substituting  (7)  for  the  partials  in  (2),  ve  find  the  coefficients  of  31/3^ 
combine  to  yield  the  le ft-band-s ide  (LHS)  of  the  expression 

[coefficient]  cosB  cos*  i*  -  sinB  l  -  ^  (l  -  )S  (8) 

[  of  31/3*  j  3r  3b  R  V  cos'*/ 


( Pape  r  1 ) 


(Paper  II) 
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The  partials  3^./3r  and  3<f//3h  are  obtained  from  ( *0  :  -sinv/R  and  cosi'/B,  respectively.  Next, 

C5)  Is  used  to  replace  cosBcosX  in  the  LHS  of  (6)  which  leaves  only  sinB  as  a  remaining  Paper  I 
coordinate.  Substituting  (5)  and  (6)  for  X  in  1  =  cos2X  ♦  sin2X  results  in  a  quadratic  ex¬ 
pression  which  can  be  solved  for  sinB, 


sinB  »  -cose  sinii 


cos* 


sine 


2  ■ 

cos^i  \ 

COS4-^! 


(9) 


The  ambiguity  ir.  sign  is  resolved  by  considering  the  special  case  of  W  0  in  (5)  and  (6)  which 
yields 

6  «  e  -  *  for  i  -  90°  .  (10) 


Relation  (10)  is  preserved  for  the  positive  root  in  (9).  The  equivalence  of  (3)  then  trivially 
follows  by  substituting  (5)  and  ( 10  with  ♦  sign)  into  its  LHS.  Alternatively,  one  can  demon¬ 
strate  the  equivalence  of  (8)  by  solving  the  spherical  triangle  in  Fig.  2  of  Paper  I. 

The  equivalence  of  the  coefficients  3I/3R  and  3I/ct  can  be  demonstrated  in  the  same  manner 
as  that  used  above  for  31/3^.  The  proof  is  less  cumbersome  than  for  31/3^  and  will  not  be  given 
here . 

As  a  final  comment,  we  note  that  although  the  equivalent  expressions  (2)  and  (3)  both  rep¬ 
resent  -he  general  inversion  of  the  brightness  integral,  the  values  of  all  three  partial  deri¬ 
vatives  cannot  be  obtained  by  a  single  observer  moving  in  solar  orbit.  Only  one  (of  the  two) 
spatial  derivatives  can  be  obtained  from  a  spacecraft.  This  limits  the  practical  use  of  the 
inversion  formula  to  "loos"  directions  confined  to  a  single  plane,  as  discussed  in  Paper  I  and 
by  Dumont  et  al .  (1979),  rather  than  the  full  steradians. 
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ABSTRACT 

In  an  earlier  paper  Sparrow  ei  al.  (1976)  found  the  polarized  bright¬ 
ness  of  zodiacal  light  to  have  solar  color  at  five  sky  positions  for  which 
there  were  fixed-position  observations  from  Skylab:  north  celestial  pole, 
south  ecliptic  pole,  vernal  equinox,  and  two  places  near  the  nor.th  galac¬ 
tic  pole.  The  brightness  and  degree  of  polarization  of  zodiacal  light  at 
these  sky  positions  are  derived  using  Pioneer  10  observations  of  back¬ 
ground  starlight  from  beyond  the  asteroid  belt  (Weinberg  et  al.,  1974; 
Schuerman  et  al.,  1976)  and  the  assumption  that  the  zodj  cal  light  is 
also  solar  color  in  total  light. 

Ten-color  observations  of  sky  brightness  and  polarization  were  made 
of  portions  of  the  antisolar  hemisphere  from  Skylab  (Weinberg  et  al., 

1975;  Sparrow  et  al. ,  1977).  Three  independent  quantities  were  determined 
directly  from  the  measurements:  total  brightness  (Be),  polarized  bright¬ 
ness  (Bp),  and  orientation  of  the  plane  of  polarization  (y) ,  the  last  more 
precisely  referred  to  as  azimuth  of  vibration.  These  three  quantities  are 
related  to  the  Stokes  parameters  I,  Q,  and  U,  respectively. 

Except  in  regions  relatively  close  to  the  Milky  Way,  the  observed 
direction  and  amount  of  polarization  (y,Bp)  in  these  space  measurements 
can  generally  be  attributed  entirely  to  zodiacal  light.  To  isolate  the 
zodiacal  light  brightness,  it  is  necessary  to  remove  the  brightness  con¬ 
tribution  of  all  discrete  or  "resolved"  stars  in  each  6°  diameter  field 
of  view  (FOV)  and  to  subtract  the  background  starlight  (integrated  star¬ 
light,  diffuse  galactic  light,  light  from  extragalactic  sources).  Thus, 
the  zodiacal  light  total  brightness  (Bz)  is  obtained  as  the  residual  after 
a  several-step  subtraction  process.  This  subtraction  process  has  been  the 
largest  single  source  of  error  in  determining  the  brightness  of  zodiacal 
light  in  observations  from  space  and,  even  more  so,  in  observations  from 
the  ground. 

Skylab  observations  included  sky-scanning  and  fixed-position  measure¬ 
ments,  the  latter  at  the  north  celestial  pole,  south  ecliptic  pole,  and 
vernal  equinox,  and  at  two  positions  near  the  north  galactic  pole.  A 
number  of  sequences  of  filters  (Table  1)  were  used  to  observe  these 
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positions,  one  sequence  taking  2  minutes.  In  the  following,  we  briefly 
outline  our  method  for  evaluating  discrete  and,  especially,  background 
starlight  at  the  Skvlab  wavelengths,  so  as  to  derive  the  total  brightness 
of  zodiacal  light  at  these  sky  positions.  The  method  will  be  developed 
more  fully  elsewhere,  in  conjunction  with  its  use  to  isolate  the  zodiacal 
light  in  Skylab  sky-scanning  observations. 

Omitting  absolute  calibration  and  the  conversion  of  analog  voltages 
to  relative  brightness  and  polarization,  reduction  of  the  Skylab  obser¬ 
vations  involves  the  following  major  steps: 

1.  Determination  of  instrument  pointing.  Pointing  directions  are  ob¬ 
tained  by  combining  vehicle  attitude  data  with  telemetry  data  on  position 
in  the  instrument's  reference  system  and  using  data  derived  from  photo¬ 
graphs  with  a  bore-sighted,  16mm  camera.  Where  there  are  frames  contain¬ 
ing  at  least  three  or  four  identifiable  stars,  pointing  is  generally  known 
to  at  least  ±0.2° . 

2.  Determination  of  the  brightness  contribution  of  stars  resolved  by  the 
instrument.  A  combined  star  catalog  was  generated  by  editing  and  merging 
data  in  the  U.S.  Naval  Observatory  Photoelectric  Catalog  (Blanco  et  al .  , 
1968),  the  Yale  Catalogue  of  Bright  Stars  (Hoffleit,  1964),  and  the  Moscow 
General  Catalog  of  Variable  Stars.  The  Skylab  spectral  transmittance 
functions  were  convolved  with  40  stellar  spectral  energy  distributions 
selected  from  the  five  luminosity  classes  given  in  Mitchell  and  Johnson 
(1969).  These  were  used  to  determine  the  relationship  of  the  instrument 
magnitude  to  the  Johnson  V  magnitude  as  a  function  of  the  B-V  color  index 
and  spectral  luminosity  class.  This  relationship  was  then  used  to  convert 
the  Johnson  V  magnitudes  of  the  combined  star  c,  log  to  instrument  mag¬ 
nitudes  for  each  Skylab  filter. 

3.  Separation  of  background  starlight  and  zodiacal  light  at  Skylab  wave¬ 
lengths.  The  Pioneer  10/11  Imaging  Photopolarimeters  were  used  to  period¬ 
ically  measure  skv  brightness  and  polarization  in  the  blue  (B/W  3950A- 
4850&)  and  red  (5900A-6900&)  at  heliocentric  distances  beyond  1.002  AU . 

A  two-color  map  of  background  starlight  over  the  sky  has  just  been  com¬ 
pleted  using  Pioneer  10  observations  from  beyond  the  asteroid  belt 
(Weinberg  et  at.,  1974),  where  the  zodiacal  light  was  found  o  be  negli¬ 
gible  compared  to  the  background  starlight  (Hanner  et  al.,  1974;  Schuerman 
et  at.,  1977).  To  create  this  map,  the  contributions  of  resolved  stars 
(stars  brighter  than  V  magnitude  =  6.5)  were  removed  from  each  of  60,000 
FOV ' s  before  merging  data  from  sky-mapping  observations  made  at  six  helio¬ 
centric  distances  between  3.27  AU  and  5.15  AU.  These  data  are  used  to 
separate  background  starlight  and  zodiacal  light  in  those  Pioneer  10/11 
observations  which  contain  both;  i.e.,  in  observations  at  heliocentric 
distances  less  than  approximately  3  AU.  The  Pioneer  background  starlight 
data  can  also  be  used  directly  in  any  other  observations  having  similar 
spatial  and  discrete  star  resolution  -  and  similar  wavelengths  (blue  and 
red)  . 


Table  1  illustrates  our  use  of  Pioneer  red  background  starlight 
(effective  wavelength  6420A)  directly  with  data  at  Skylab  wavelength 
6427A.  Also  shown  is  the  method  used  for  other  wavelengths  in  Skylab 
fixed-position  observations  at  the  north  celestial  pole. 
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Table  1.  Observed  and 

Der ived 

Results 

at  the 

North  Celestial 

Pole 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Filter  central 

wavelength 

X 

Bn 

Bt 

Bd  s 

B t —  Bd  s 

Bbs 

BZ 

400LA 

24!8 

23.8 

192 

28.0 

164 

[87] 

77 

4  /  48A 

27.8 

18.5 

154 

14.3 

140 

63 

77 

5068A 

28.0 

19.4 

164 

13.7 

150 

73 

77 

5294A 

27.7 

18.3 

163 

12.5 

151 

74 

77 

5562A 

27.3 

19.3 

179 

11.2 

168 

[91] 

77 

606  3A 

26.7 

17.9 

164 

11.3 

153 

76 

77 

6286X 

27.0 

19.2 

221 

11.9 

209 

[132] 

77 

6  4  2  7  A  . 

26.6- 

*  • 

18.7  • 

•  •  167  •  •  • 

11.6-  • 

.  .  155 - 

•  78-  •  • 

•  77 

7093& 

28.4 

18.9 

172 

11.3 

161 

84 

77 

8160& 

29.4 

15.4 

187 

11.5 

176 

[99] 

77 

Mean  Values 

27  ?4± . 

,  1 

18. 8±. 

,2 

77 

Each  value  of  y, 

Bp ,  and 

Bt  is 

the  mean 

of  up 

to  15  individual 

measure 

ments . 

Table  1  columns  and  discussion: 

(2)  Measured  orientations  of  the  plane  of  polarization;  these  values 
correspond  closely  to  the  electric  vector's  being  perpendicular  to 
the  scattering  plane. 

(3)  Polarized  brightness  (Bp),  in  Sjq(v)  units.  The  Sjq(V)  unit  is  re¬ 
ferred  to  the  sun;  the  brightnesses  in  this  unit  will  have  the  same 
value  at  all  wavelengths  for  radiation  of  solar  color.  Absolute  cal¬ 
ibrations  at  400lX  and  8160&  appear  to  be  high  and  low,  respectively, 
as  noted  in  a  first  analysis  of  these  data  by  Sparrow  et  at.  (1976). 
The  mean  value  is  determined  using  8  colors  only. 

(4)  Total  observed  brightnesses  (Bt) .  These  include  discrete  and  back¬ 
ground  starlight,  zodiacal  light,  and  airglow  emission;  the  last  at 
5562A,  6286A,  and  probably  8160A. 

(5)  Discrete  starlight  (Bds) . 

(6)  Total  brightness  (Bt)  minus  discrete  starlight  (Bds)-  Except  for  the 
aforementioned  airglow,  these  differences  correspond  to  background 
starlight  plus  zodiacal  light. 

(7)  and  (8):  Background  starlight  (B^s)  and  zodiacal  light  (Bz) ,  respec¬ 

tively.  For  Bbs  Skylab  642 7 A  is  equivalent  to  Pioneer  red,  which  we 
found  (column  (7))  to  be  78  Sjo(V).  Thus,  Bt-Btjs-B^)S=BZ=77  Siq(V). 

As  seen  in  column  (3)  (see,  also.  Sparrow  et  at.,  1976),  Bp  is  close 
to  solar  color.  If  we  assume  that  Bz  is  also  solar  color,  then  Bz= 

77  Sjo(V)  aC  all  Skylab  wavelengths  (column  (8)).  B^g  at  the  other 
Skylab  wavelengths  is  obtained  from  the  corresponding  differences, 
(6)-(8).  Due  to  calibration  difficulties  and  airglow  emission,  the 
bracketed  values  of  B^g  are  uncertain.  From  Pioneer  blue  data 
(effective  wavelength  4360A)  we  find  B^,s=57  Sjq(V),  in  good  agreement 
with  the  trend  in  Skylab-der ived  values  of  B^s;  i.e.,  the  results  are 
consistent  with  the  assumption  of  solar  color  for  Bz. 

Table  2  gives  the  results  obtained  by  applying  this  method  to  all 
fixed-position  observations.  Each  value  of  y  and  Bp  is  a  mean  value  ob¬ 
tained  as  shown  in  Table  1.  Only  part  of  the  data  was  available  for  the 
earlier  analysis  of  the  Bp  data  by  Sparrow  et  al .  (1976),  which  also  found 
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Table  2.  Best  Estimate  of  Mean  Zodiacal  Light  -  Fixed  Position 
Observations,  Skylab  Mission  SL-2 


Coordinates 

Zodiac 

a  1 

Light 

Target 

Equal 

torial 

Ecliptic 

a 

6 

e 

6 

X 

bp 

Bz 

p 

NCP 

174?2 

85“8 

67?2 

65 

“8 

27! 

4±. 

1  18. 

8± 

.2 

77 

.244 

SEP 

76.2 

-66.6 

89.8 

-84 

.5 

30. 

8± . 

3  11. 

4t 

.5 

50 

.228 

VE 

350.6 

1.0 

89.9 

4 

.6 

29. 

6± . 

2  25. 

6± 

.4 

168 

.152 

NCP1 

190.9 

27.6 

95.4 

29 

.4 

29. 

8i . 

4  17. 

Oi 

110 

.154 

NGP2 

189  S 

27.4 

92.4 

28 

.8 

29. 

51. 

2  17. 

2  ± 

.2 

112 

.154 

Bn  to  be 

solar  color.  Availabil 

itv 

cf 

all 

of 

the  da 

ta 

f  or 

our 

analysi s 

and  our 

use  of  a 

different  solar 

-based 

S1  0 

(V) 

unit  a 

ccount 

for 

the  sligh 

lower  values  of  Bp  found  here.  In  this  paper  and  in  our  studies  of  Pioneer 
10/11  data,  we  follow  the  Sjq(V)  recommendations  of  Sparrow  and  Weinberg 
(1976).  The  degree  of  polarization,  p,  is  derived,  not  measured.  It  is 
obtained  from  Bp/Bz.  The  measured  and  assumed  solar  colors  for  Bp  and  Bz, 
respectively,  mandate  that  p  be  grey  or  independent  of  wavelength.  The  p 
data  reconfirm  our  earlier  result  (Weinberg  and  Hahn,  1976)  and  that  of 
Dumont  and  Sanchez  (1976)  that  p  is  greater  off  the  ecliptic  than  on  the 
ecliptic  at  these  elongations.  As  noted  earlier,  this  methodology  is  also 
being  used  in  the  reduction  ol  the  Skylab  scanning  data. 

We  thank  Dr.  Frank  Giovane  for  developing  the  Skylab  star  subtraction 
methodology  and  Susan  Darbys'nire  for  assisting  in  the  reduction  of  these 
data.  This  work  is  supported  by  NASA  grant  NSG  8040. 
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ABSTRACT 

The  classical  restricted  three-body  problem  is  generalized  to  in¬ 
clude  the  force  of  radiation  pressure  and  the  Poynting-Rober tson  effect, 
The  positions  of  the  Lagrangian  points  L^  and  L5  are  found  as  functions 
of  6,  the  ratio  of  radiational  to  gravitational  forces.  The  Poynting- 
Robertson  effect  renders  the  L^  and  L5  points  unstable  on  a  time  scale 
(T)  long  compared  to  the  period  of  rotation  of  the  two  massive  bodies. 
For  the  solar  system,  T  is  given  by  T  =  [ ( 1— S )7°/B]  544  a2  yr  where  a 
is  the  separation  between  the  Sun  and  the  planet  in  AU.  Implications 
for  space  colonization  and  a  mechanism  for  producing  asymmetries  in  the 
interplanetary  dust  complex  are  discussed;  testing  the  latter  may  be 
possible  from  the  Zodiacal  Light/Background  Starlight  Experiment  aboard 
the  International  Solar  Polar  Mission  spacecraft  to  be  launched  in  1983 


A  steady  state  solution  of  the  motion  of  a  B  particle  (a  particle 
whose  surface  to  volume  is  such  that  the  force  of  radiation  pressure  is 
non-negligible)  in  the  gravitational  field  of  two  larger,  orbiting 
bodies  of  masses  M-^  and  M2  is  sought.  It  is  assumed  that  M^  and  M2  are 
in  circular  orbits  about  their  common  center  of  mass,  that  their  angular 
frequency  (Cl)  is  given  by  fl2  =  G(M|  +  M2)/a3,  an d  that  the  6  particle 
remains  in  the  plane  in  which  M^  and  M2  revolve.  In  the  rotating  refer¬ 
ence  frame  whose  origin  is  the  center  of  mass  and  whose  angular  frequency 
is  Cl,  the  two  components  of  the  equation  of  motion  of  a  B  particle  with 
coordinates  (X, Y)  are 


X  2Y  ..  (1-Bi )u (X+u-1)  (1-8?  ) (1-u) (X+u ) 

_2  _  _  =  x - ^ 


Y  .  2ft  r,  (1-Bi )u  ( 1-8- )  (l-u)i„ 
T2  +  q  =  l1  H  r  q  >Y 


(1) 

(2) 
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M,  is  located  at  X  =  1-u ,  M2  is  located  at  X  =  -p  and  u  =  M^/(M^+M7). 

The  distances  from  the  point  (X,Y)  to  the  centers  of  mass  M-^  and  M2  are 
rj  =  |(x+u-l)2  +  Y2 } • 2  and  r2  =  j(x+p)2  +  Y2|’2.  respectively.  The  two 
terms  on  the  left-hand  side  of  equations  (1)  and  (2)  represent  the 
acceleration  and  the  coriolis  force  experienced  by  the  £  particle;  the 
centrifugal  acceleration  has  been  placed  on  the  right-hand  side  (RHS) . 

The  last  two  terms  on  the  RHS ' s  of  equations  (1)  and  (2)  represent  the 
components  of  the  "effective"  gravitational  forces.  These  terms  differ 
from  those  of  the  classical  problem  in  that  each  component  has  been  re- 
auaca  bu  a  factor  of  the  fern  1-6. 

The  equilibrium  points  L4  and  L5  are  found  by  setting  the  RHS's  of 
equations  (1)  and  (2)  equal  to  zero  with  Y  /  0.  The  solution  is 

r  1  =  ( 1— S 1  )’3  =  6 1  ,  r2  =  (1-62)  3  =  &2  •  (3) 

From  solution  (3)  it  follows  that  the  L4  and  L5  points  exist  only  if  all 
of  the  following  conditions  are  satisfied:  6i+<$2  >  1,  0  <  <  1,  and 

0  <  {j  <  1.  These  results  are  plotted  in  Fig.  1  which  shows  the  positions 
of  the  and  L5  points  as  a  function  of  gj  and  62  (and  equivalently, 

S [  and  62)-  For  application  to  the  solar  system,  62  vanishes  and  62  is 
unity;  the  locus  of  these  solutions,  parameterized  by  gj ,  is  the  heavy 
arc  shown  in  Fig.  1.  (Henceforth,  the  subscripts  on  Bj  and  62  will  be 
omitted).  Furthermore,  these  solutions,  when  applied  to  the  solar 
system  (u-1)  are  stable  in  the  sense  that  a  small  displacement  of  the 
3  particle  from  equilibrium  will  result  in  oscillations  whose  amplitude 
does  not  increase  with  time. 


FIG.  1.  The  possible  positions  of 
the  equilibrium  points.  The  heavy 
line  applies  to  systems  with  only 
one  luminous  body,  M, . 


The  results  mentioned  thus 
far  have  an  interesting  appli¬ 
cation  for  artificial  satellites 
and  future  space  colonization. 

It  has  been  suggested  ( e.g .  , 
O'Neill,  1974)  that  the  classical 
triangular  points  of  the  Sun- 
Jupiter  or  Sun-Earth  system 
would  be  a  convenient  site  to 
locate  future  space  colonies. 
Present  day  technology  permits 
the  construction  of  spacecraft 
whose  value  of  6  can  be  made  or 
varied  in  the  range  0<g<l.  A 
space  colony  with  a  solar  facing 
adjustable  "sail"  could  thus 
"park"  (i.e.,  remain  stationary 
with  respect  to  the  Sun  and  a 
planet  without  expending  large 
amounts  of  energy)  not  only  at 
the  classical  triangular  positions 
but  at  any  heliocentric  distance 
out  to  that  of  the  planet. 
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When  terms  representing  the  Povnting-Robertson  effect, 


-Bnafi 


P  = 
v  cr 


H 


1  L 
-Buafi  T  Y 


x  +  X+y-1  +  p-  (  -(X+y-1)  + 


(4) 

(5) 


are  included  on  the  RHS ' s  of  equations  (1)  and  (2),  the  problem  can  still 
be  solved  analytically  (to  first  order  terms  in  afi/c)  in  the  same  manner 
as  before.  The  locations  of  the  equilibrium  points  differ  in  this 
treatment  only  by  small  corrections  of  the  order  of  afi/c.  However,  all 
solutions  are  unstable  when  the  Povnting-Robertson  effect  is  taken  into 
account.  The  time  scale  of  this  instability,  when  applied  to  the  solar 
system,  is  given  by 


T 


(l-g)?-ica; 


3BM  G 
o 


(i-e)?3 

6 


544  a" 


yr 


(a  in  AU ' s) . 


(6) 


A  convenient  and  representative  value  for  g  of  0.57  makes  (1 -6)JVB  =  1. 
Thus,  for  example,  the  Sun-Jupiter  system  (a  =  5.2  AU)  results  in 
T  =  14,000  yr.  This  means  that  if  a  particle  with  g  =  0.57  is  displaced 
slightly  from  its  equilibrium  position,  it  will  oscillate  about  that 
position  and  the  amplitude  of  that  motion  will  increase  by  a  factor  of 
e  in  14,000  yr.  (A  more  complete  derivation  of  these  results  is  given 
by  Schuertnan,  1980). 


It  is  not  at  all  evident  that  the  above  mechanism  could  lead  to  a 
significant  increase  in  the  density  of  interplanetary  g  particles  along 
the  arc  schematically  shown  as  a  heavy  line  in  Fig.  1.  Other  forces  act 
on  these  6  particles  which  cannot  be  treated  by  the  analytical  methods 
employed  here.  Interactions  with  the  gravitational  fields  of  intervening 
planets  would  leave  large  gaps  in  the  arc  of  8  particles.  Jupiter,  because 
of  its  great  mass,  would  be  particularly  effective  in  disrupting  the  dust 
in  arcs  associated  with  other  planets;  conversely,  dust  in  the  Sun-Jupiter 
arc  would  be  the  least  affected  by  planetary  encounters.  The  Lorentz 
force  due  to  the  solar  magnetic  field  also  acts  on  (charged)  6  particles. 
This  conservative  force  is  roughly  40%  that  of  the  solar  gravitation 
(Greenberg  and  Schuerman,  1978).  The  uncertainties  in  the  charges  of  the 
6  particles  and  in  the  long  term  fluctuations  and  reversals  of  the  solar 
magnetic  field  at  distances  greater  than  1  AU  make  the  inclusion  of  the 
Lorentz  force  in  this  problem  particularly  difficult.  Finally,  the  ions 
in  the  solar  wind  impinging  on  the  6  particles  and  the  Coulomb  interaction 
between  solar  wind  electrons  and  the  S  particles  contribute  non-conser¬ 
vative  drag  forces  which  together  are  estimated  to  be  about  30%  of  the 
Poynting-Robertson  effect  (Bandtrraann,  1967;  Misconi,  1976a).  These 
forces  must  surely  operate  to  shorten  the  time  scale  (T)  in  equation  (6). 

Even  if  planetary  encounters  do  not  disrupt  the  entire  arc,  and  even 
if  the  large  but  conservative  Lorentz  force  averages  out  to  a  small  effect 
over  many  solar  cycles,  there  still  remains  the  need  for  a  replenishment 
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mechanism  if  an  increase  in  the  density  of  6  particles  along  a  planetary 
arc  is  to  be  expected  on  a  "steady  state"  basis.  The  major  source  of 
replenishment  for  the  interp ianei_ary  uust  complex  is  thought  to  be 
cometarv  debris  (Whipple,  1955  and  1967;  Dohnanyi,  1970,  1972,  and  1973). 
Comets  shed  both  6  particles  (some  with  6<1)  and  larger  (il  pm)  particles; 
light  scattering  by  the  latter  is  the  main  contribution  to  the  zodiacal 
light.  These  larger  particles  ultimately  fragment  through  self-collisions 
(Zook,  1975)  o~  possibly  by  rotational  bursting  (Paddack  and  Rhee,  1975; 
Misconi,  1976b).  It  is  not  evident  whether  or  not  these  sources  of  6 
particles  provide  a  wide  enough  distribution  in  real  space  and  velocity 
space  60  as  to  deposit  a  significant  fraction  of  "zero-velocity"  S 
particles  at  their  appropriate  equilibrium  positions  in  times  short 
compared  to  T  in  equation  (6). 

Fortunately,  measurements  made  during  the  International  Solar  Polar 
Mission  may  determine  whether  or  not  this  mechanism  produces  any  "dust 
arcs"  in  the  solar  dust  complex.  A  zodiacal  light  photopolarimeter 
(jointly  designed  by  teams  at  the  Ruhr-Univer sity ,  Bochum,  FRG  and  the 
State  University  of  New  York  at  Albany)  will  map  the  brightness  and 
polarization  of  the  zodiacal  light  throughout  the  four  and  one-half  year 
mission  to  Jupiter  and  back  over  the  poles  of  the  Sun.  During  the  two 
periods  of  solar  polar  passage,  the  spacecraft  will  be  in  the  center  of 
the  solar  system  but  at  1.7  AU  above  (and  below)  the  ecliptic  plane. 

From  these  vantage  points,  isophotes  of  the  zodiacal  light  may  reveal  any 
enhancements  in  the  spatial  density  of  the  dust  along  arcs  associated 
with  planets. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  (Contract  F49620-78-C-0013)  and  by  the  National  Aeronautics 
and  Space  Administration  (JPL  Contract  955137). 
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ABSTRACT 

Traditionally,  earth-based  observations  of  the  zodiacal  light  (ZL) 
require  two  assumptions  for  further  analysis:  (Al)  the  dust  density  (n) 
is  a  power  of  heliocentric  distance  (R) ,  n  <*  R-v;  (A2)  the  nature 
(scattering  cross  section,  a)  of  the  dust  is  independent  of  location, 
a (r  ,h , 6 ) =a (e ) .  Observations  from  Pioneer  10  do  net  verify  these  assump¬ 
tions. 


Both  two-dimensional  inversion  techniques  (Schuerman,  1979a)  and 
more  traditional  methods  of  analysis  (Hanner  ei  al.,  1976)  have  been 
performed  on  portions  of  the  ZL  observations  made  from  Pioneer  10/11 
spacecraft.  In  developing  the  three-dimensional  form  of  the  ZL  inversion 
(Schuerman,  1979b;  henceforth  called  Paper  I),  it  was  realized  that  this 
more  general  technique  should  be  applied  to  the  Pioneer  data  because  the 
velocity  vectors  of  the  spacecraft  were  not  confined  to  the  ecliptic 
plane.  The  previously  mentioned  analyses  examined  data  with  viewing 


directions  parallel  to  the  ecliptic;  the 


r,h  plane  and  contains  S ,  the  direction 
of  motion  of  the  spacecraft. 

71 


three-dimensional  inversion 

requires  viewing  along  a  great 
circle  (see  Fig.  1)  containing 
(1)  the  line  (n)  which  is  both 
normal  to  the  observer-Sun 
line  and  parallel  to  the 
ecliptic  and  (2)  the  space¬ 
craft  velocity  vector  (S) . 

This  great  circle  is  called 
the  inversion  plane  in  Paper 
I  where  its  geometric  con¬ 
struction  and  the  notation 
used  here  are  defined  in  more 
detail.  The  angle  $  is  the 
inclination  of  the  inversion 
plane  to  the  ecliptic,  and 
viewing  in  the  direction  of  n 
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Table  1 


COORDINATES  OF  SPACECRAFT  POSITION  (r,h  in  AU)  AND  VALUES  OF 
THE  INCLINATION  OF  THE  INVERSION  PLANE  ($  in  deg)  AS  A 
FUNCTION  OF  HELIOCENTRIC  DISTANCE  (R  in  AU) 


R: 

1 . 010 | 1 . 150 

1.171 

1.349 

1.453 

1.652 

1.865 

2.294 

2.413 

2.467 

2.641 

2.938 

r: 

1.010  1.150 

1.171 

1.349 

1.452 

1.651 

1.864 

2.292 

2.412 

2.465 

2.640 

2.937 

h: 

-.0091-. 030 

-.032 

-.045 

-.051 

-.067 

-.079 

-.081 

-.082 

-.085 

-.089 

4> : 

— 12.9  j— 8.60 

-8.04 

-4.11 

-3.29 

HBWE1 

BSB 

-1.01 

-0.82 

-0.51 

corresponds  to  an  elongation  of  e  =  90°.  The  observer  (spacecraft)  is 
located  by  cylindrical  ecliptic  coordinates  (r,h).  Table  1  lists  the 
values  of  $  and  h  as  a  function  of  R(=r)  for  Pioneer  10.  Notice  that 
$  is  small  as  is  the  ratio  h/r.  Nevertheless,  we  have  isolated  those 
Pioneer  data  having  views  along  the  inversion  plane  and  have  analyzed 
those  data  via  both  inversion  and  traditional  methods.  The  former  anal¬ 
ysis  will  be  published  elsewhere.  The  blue  channel  brightness  data  in 
units  of  Siq(V)  are  shown  in  Figure  2.  The  total  ZL  brightness  (Z)  as 
a  function  of  e  was  derived  for  12  distinct  heliocentric  distances  in 
order  to  generate  this  Z(R,e)  surface.  The  values  of  Z  have  been  smoothed 
somewhat  by  a  curve-fitting  procedure  which  also  yields  a  measure  of  the 
statistical  uncertainty  (o2)  in  Z  at  each  point  on  the  surface.  Now 
we  can  apply  the  traditional  method  of  ZL  analysis:  for  each  value  of 
e,  we  fit  a  power  law  of  the  form  Z  «  Ra  where  a  =  -(1+v) ;  see,  for 
example,  Leinert,  (1975).  If  assumptions  A1  and  A2  are  correct ,  a  single 
value  of  a  should  apply  to  all  values  of  elongation. 

The  logZ  vs  logR  plot  in 
Figure  3  shows  a  typical  example 
of  the  fitting  procedure  used  to 
determine  a,  the  slop*  of  the 
assumed  straight  line,  for 
e  =  115°.  Each  data  point  on 
the  fit  is  represented  by  two 
numbers  of  the  form  Z/oz.  The 
vertical  error  bars  represent 
the  corresponding  uncertainty 
in  logZ ;  o^ogz  =  ln(10)  oz/Z. 
Notice  that  even  though  some 
data  points  have  remarkably  well 
determined  (statistically  speak¬ 
ing)  values  (3/2,  4/1,  etc...), 
the  innate  mathematical  proper¬ 
ties  of  the  assumed  Z  a  ra  re¬ 
lation  produce  large  uncertain¬ 
ties  in  logZ  .when  Z  is  small. 
(Also,  Z  itself  becomes  more 
difficult  to  separate  from  the 
background  starlight  for 


Fig.  2.  The  brightness  (Z)  of  the 
ZL  as  a  function  of  heliocentric 
distance  (R)  and  elongation  (e). 
Viewing  directions  are  confined 
to  the  inversion  plane. 
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Fig.  3.  LogZ  vs  LogR  for  e  =  115°.  The  slope  of  the  fitted  line  is 
a  =  -(1+v).  The  number  pairs  represent  the  data  in  the  form  Z/o2.  The 
vertical  bars  are  the  uncertainty  in  logZ  given  by  ln(10)  o2/Z. 

observations  made  at  large  c  and  R.  In  this  regard,  ZL  studies  made 
from  outgoing  spaceprobes  require  very  different  data  handling  tech¬ 
niques  than  those  used  for  solar  approaching  probes  like  Helios.') 
Following  standard  statistical  practice,  we  have  weighted  each  data 
point  by  ogz)-^  in  determining  the  slope  a.  For  the  example  shown 
in  Figure  3,  a  =  -3.6  ±  1.4. 

Fitted  values  of  a  were  obtained  in  the  above  manner  for  all 
70°<  e  <175°  in  intervals  of  hz  =  5°.  Figure  4  depicts  the  result  of 
those  fits  as  a  function  of  e.  The  vertical  bars  indicate  the  uncertain¬ 
ty  in  a.  Earlier  and  independent  Pioneer  results  by  Hanner  et  at.  (1976) 
are  shown  as  the  shaded  rectangle.  Their  estimates  of  a  are  in  agree¬ 
ment  with  those  found  here  over  the  limited  range  of  elongation  that 
they  considered.  However,  we  find  a  systematic  decrease  in  a  for 
85°$  c  Sl25°.  It  is  significant  that  near  e  =  90°  the  "dip"  in  a  can 
not  be  attributed  to  viewing  in  the  inversion  plane;  the  line  of  sight 
is  exactly  parallel  to  the  ecliptic  for  e  =  90°.  Yet,  it  is  near  this 
direction  that  our  results  differ  most  from  the  a  =  -2.3  value  reported 
from  Helios  (Link  et  at.,  1976).  As  a  final  check  on  our  results,  we 
similarly  analyzed  the  in-ecliptic  Pioneer  data  previously  published 
by  Schuerman  (1979a) .  The  paucity  of  that  data  restricts  the  analysis 
to  e >1 10 ° .  Nevertheless,  the  same  general  trend  was  found  as  denoted 
by  the  x's  in  Figure  4.  The  systematic  drop  of  a  towards  z  *=  100°  is 
evident,  and  the  scatter  in  those  results  is  well  contained  within  the 
envelope  of  our  error  bars. 

Although  we  are  pursuing  other  explanations,  some  of  which  address 
the  fact  that  h/R  is  not  exactly  zero,  we  tentatively  maintain  that 
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Fig.  4.  Values  of  the  index  a  (in  Z  “  Ra)  as  a  function  of  elongation 
(e).  The  circles  and  associated  error  bars  refer  to  the  data  shown  in 
Fig.  2.  Data  points  marked  by  x's  refer  to  in-ecliptic  Pioneer  data 
previously  presented  by  Schuerman  (1979a).  The  shaded  rectangle  approxi¬ 
mates  the  earlier  and  independent  Pioneer  results  of  Hanner  et  al. 

(1976).  If  the  scattering  properties  of  the  dust  are  independent  of 
position,  a  should  have  the  same  value  for  all  s. 


observational  evidence  now  exists  that  either  assumption  Al  or  A2  or 
both  are  incorrect.  In  fact,  assumption  A2  has  no  theoretical  support 
whatsoever.  On  the  contrary,  most  forces  tend  to  separate  the  dust  part 
icles  (as  a  function  of  R)  according  to  their  charge  to  mass  ratios, 
surface  area  to  volume,  etc.  If  the  results  presented  here  are  verified 
it  means  that  optical  observations  from  space  can  truly  begin  to  address 
the  dynamics  governing  the  solar  system  dust  complex. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  contract  F49620-78-C-0013. 
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ABSTRACT 

Analysis  of  zodiacal  light  observations  from  Mt.  Haleakala,  Hawaii 
show  that  the  symmetry  plane  of  the  zodiacal  cloud  near  1  A.U.  is  close 
to  the  invariable  plane  of  the  solar  system.  Since  the  symmetry  plane 
of  the  inner  zodiacal  cloud  is  close  to  the  orbital  plane  of  Venus 
(Misconi  and  Weinberg,  1978;  Leinert  et  al. ,  1979),  we  suggest  that  the 
symmetry  plane  changes  inclination  with  heliocentric  distance. 

INTRODUCTION 

The  location  in  space  of  the  plane  of  maximum  dust  density  (symmetry 
plane)  of  the  zodiacal  cloud  has  important  implications  for  the  dynamics 
of  solid  particles  in  the  solar  system.  These  implications  are  related 
to  questions  such  as  the  gravitational  perturbations  of  the  planets  on 
the  dust,  the  effects  of  solar  magnetic  forces,  and  the  general  dynamical 
history  of  dust  ejected  from  comets.  Helios  1  and  2  zodiacal  light 
measurements  near  the  region  of  1  A.U.  (Leinert  et  al.,  1979)  find  that 
the  symmetry  plane  in  this  region  has  an  inclination  i  =  3  ±0.3°  and 
ascending  node  tJ  =  87  ±  4° ,  which  is  different  from  that  of  the  invariabl 
plane  of  the  solar  system  ( i  =  1.6°,  0.  =  107°).  At  the  same  time,  ground 
based  observations  at  Tenerife  in  the  Canary  Islands  (Dumont  and  Sanchez, 
1968)  and  satellite  observations  from  the  satellite  D2A  (Dumont  and 
Levasseur-Regourd ,  1978)  find  the  symmetry  plane  close  to  the  invariable 
plane  near  1  A.U. 

In  the  present  study  an  attempt  is  made  to  obtain  additional  data 
on  the  location  of  the  symmetry  plane  near  the  Earth’s  orbit  using 
ground-based  observations  of  the  zodiacal  light  from  Mt.  Haleakala, 

Hawaii  taken  by  Weinberg  and  Mann  (1967) . 

THE  OBSERVATIONS 

Two  nights  of  observations  were  found  to  be  suitable  for  this 
purpose:  the  evenings  of  23  and  24  February,  1968.  These  photoelectric 
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observations  were  taken  using  a  narrow  band  filter  centered  at  5300  X 
(Weinberg  and  Mann,  1967).  The  photometer  made  almucantar  scans  of  the 
evening  zodiacal  light  cone  in  5°  steps  of  elevation  from  10  to  45°. 


Fig.  1  shows  an  isophote  map  of  the  zodiacal  light  brightness  in 
ecliptic  coordinates  S)  on  February  24,  1968.  Details  of  the  data 

reduction  and  the  method  of  finding  the  photometric  axis  (dashed  line) 
were  given  earlier  (Weinberg,  1964;  Misconi,  1976  and  1978).  The  dis¬ 
placement  in  ecliptic  latitude  of  the  photometric  axis  is  found  to  be 
between  1  and  1.25°  from  the  ecliptic  plane.  In  Fig.  2,  we  show  the 
displacements  from  the  ecliptic  plane  of  the  observed  plane  (dashed  line) 
of  maximum  brightness  on  two  nights  in  February,  1968,  at  elongation 
angles  50-90°  from  the  Sun.  The  predicted  displacements  of  the  invariable 
plane  and  the  orbital  plane  of  Venus,  and  the  extension  of  the  svmnetrv 
plane  found  by  Misconi  and  Weinberg  (1978)  correspond  to  what  would  be 
observed  if  the  dust  that  contributes  most  to  the  maximum  brightness 
were  in  these  planes  (Misconi,  1977).  It  is  clear  that  the  observed 
"plane"  departs  from  the  symmetry  plane  of  the  inner  zodiacal  cloud  and 
approaches  trie  invariable  plane  near  the  Earth's  orbit:  we  find  that  the 
inclination  decreases  from  i  ~  2.7°  at  c  =  50°  to  i  -  1.8  to  2°  at  c  ~  90° 


ecliptic 


5  ic 

latitude 


Fig.  1.  Brightness  contours 
the  evening  zodiacal  light  on 
February  24,  1968  in  ecliptic 
coordinates  (>->0,£). 
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Fig.  2.  Predicted  displacements  in 
ecliptic  latitude  (solid  lines)  of 
the  orbital  plane  of  Venus,  extension 
of  the  symmetry  plane  of  the  inner 
zodiacal  light  (e  32-50°)  found  by 
Misconi  and  Weinberg  (1978),  and  the 
invariable  plane,  as  a  function  of 
elongation  angle.  The  observed 
symmetry  plane  in  this  region  is 
shown  by  the  dashed  line. 
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At  the  same  tine,  Leinert  ci  al.  (1979)  found  that  a  single  inclination, 
7=3*  0.3°,  applies  over  this  entire  range  of  elongation.  As  noted 
earlier,  Dumont  and  Sanchez  (1968)  and  Dumont  and  Levasseur-Regourd  (1978) 
find  the  symmetry  plane  near  1  A.U.  to  be  close  to  the  invariable  plane. 
Fig.  3  sketches  the  relative  positions  of  the  observed  symmetry  plane 
(dashed  line)  by  Leinert  ct  al.  (1979),  our  observed  svmnetrv  "surface" 
as  we  prefer  to  call  it  (solid  line),  the  invariable  plane,  the  orbital 
plane  of  Venus,  and  the  solar  equatorial  plane. 


Fig.  3.  A  sketch  of  the  relative  inclinations  from  the  ecliptic  plane 
as  a  function  of  elongation:  for  the  solar  equatorial  plane,  the  orbital 
plane  of  Venus  and  the  invariable  plane.  Also  shown  is  the  position  of 
the  symmetry  plane  found  by  Leinert  ez  al.  (dashed  line)  and  our  combined 
results  over  this  range  of  elongation. 


DISCUSSION  OF  RESULTS 

It  is  found  that  the  symmetry  plane  of  the  zodiacal  cloud  at  elon¬ 
gation  angles  50-90°  does  not  continue  to  be  close  to  the  orbital  plane 
of  Venus  but  rather  departs  from  it  and  approaches  the  invariable  plane. 

In  view  of  this  result  one  must  re-evaluate  the  dynamical  forces  influ¬ 
encing  the  spatial  distribution  of  the  dust.  If  one  considers  the  effects 
of  the  solar  magnetic  forces  found  by  Griin  and  Morfill  and  reported  in 
this  symposium  it  would  seem  that  the  zodiacal  cloud  should  have  one 
symmetry  plane  decreasing  smoothly  in  inclination  with  increasing  helio¬ 
centric  distance.  Our  result  shows  this  tendency,  but  additional  data 
are  needed.  Since  the  inclination  of  the  inner  part  of  the  observed 
plane  is  2.7°  and  the  inclination  of  the  solar  equatorial  plane  is  7.3°, 
the  solar  magnetic  forces  have  not  yet  achieved  a  symmetry  to  the  solar 
equatorial  plane.  This  fact  raises  questions  about  the  lifetime  of  these 
effects  and  the  origin  of  the  dust  in  the  inner  region  of  the  solar  system. 

It  is  not  possible  on  an  observational  basis  alone  to  distinguish 
between  the  gravitational  effects  of  Venus  and  solar  electromagnetic 
forces  on  the  dust.  This  comes  about  from  the  fact  that  the  solar  equa¬ 
torial  plane  and  the  orbital  plane  of  Venus  have  almost  the  same  ascend¬ 
ing  nodes  (ft  =  73°  and  ft  =  77°,  respectively).  Therefore,  the  question 
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of  what  is  influencing  the  dust  distribution  is  still  unanswered.  A 
dynamical  evaluation  of  the  gravitational  influence  of  the  inner  planets 
on  the  dust  distribution  is  needed  before  one  can  discriminate  between 
these  effects. 

New  observations  of  the  zodiacal  light  are  scheduled  from  Mt. 
Haleakala  at  c  =  50-180°  to  provide  more  information  on  the  svmmetry 
plane  in  these  regions. 

We  thank  D.  Beeson  and  Tani  Engel  for  assistance.  Supported  by 
NASA  grant  NSG  7093  and  by  Air  Force  Office  of  Scientific  Research 
Contract  F  49620-78-C-0013 . 
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DISCUSSION 

Srdn:  Electromagnetic  effects  (Griin  and  Morfill,  this  volume)  perfectly 
account  for  your  observed  variation  of  the  zodiacal  light  symmetry 
plane  with  solar  distance,  i.e.,  the  increase  of  inclination  toward  the 
Sun.  In  particular,  the  observed  longitude  of  the  ascending  node  out¬ 
side  the  Earth's  orbit  supports  the  importance  of  the  solar  equator  and 
contradicts  the  gravitational  influence  of  Venus  on  dust  particles  out¬ 
side  the  Earth's  orbit. 

Mo s com :  Electromagnetic  effects  are  certainly  one  source  affecting  the 
symmetry'  plane  of  the  zodiacal  cloud.  From  these  effects  one  would 
expect  a  smooth  decrease  in  the  inclination  of  the  symmetry  plane  away 
from  the  Sun,  whereas  our  observations  show  an  irregular  decrease  in 
inclination.  We  agree  that  the  fact  that  the  symmetry  plane  is  near  the 
orbital  plane  of  Venus  does  not  necessarily  mean  that  the  dust  is 
perturbed  gravitationally  by  Venus.  But,  unless  the  lack  of  a  gravi¬ 
tational  effect  by  Venus  is  shown  through  a  treatment  using  perturbation 
techniques  -  which  to  our  knowledge  has  not  been  done  -  gravitational 
effects  along  with  electromagnetic  effects  should  continue  to  be  con¬ 
sidered  as  potential  forces  affecting  the  symmetry  plane.  Our  obser¬ 
vations  find  the  symmetry  plane  near  the  orbital  plane  of  Venus  only  at 
heliocentric  distances  between  0.5  and  0.8  AU,  not  outside  the  Earth's 
orb i t . 
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r.-'Szh:  How  does  Dr.  Tanabe's  measurement  of  the  inclination  near  the 
Gegenschein  compare  with  measurements  out  to  Earth? 

Vic  nt ;  Dr.  Tanabe’s  results  show  that  the  symmetry  plane  outside  the 
Earth's  orbit  decreases  in  inclination  to  a  value  of  0.5°  and  an 
ascending  node  close  to  that  of  the  orbital  plane  of  Venus;  i.e.,  a 
general  decrease  in  the  inclination  away  from  the  Sun.  Dr.  Tanabe  also 
suggests  that  the  combined  masses  of  the  Earth  and  Venus  are  affecting 
dust  outside  the  Earth's  orbit. 

ru!n:r :  A  calculation  has  recently  been  completed  (Burns,  J.A.,  Hamill, 
P.,  Cuzzi,  J.N.  and  Durisen,  R.:  1979,  Astron.  J.  in  press)  that  might 
pertain  to  the  question  of  gravitational  perturbation  of  the  symmetry 
plane.  We  studied  analytically  the  out-of-plane  perturbations  of  a 
particle  in  Satums's  ring  as  affected  by  an  inclined  satellite  of 
Saturn  and  the  planet's  oblateness.  We  found  that  the  ring  particles 
suffered  short  -  period  oscillations  of  amplitude  (mpa r t /M) ( a/apar r )  " 
a(sin  ipart)  as  the  particle  passed  under  the  perturbing  satellite.  The 
long-  period  perturbations  were  smaller  than  this  but  of  the  same  form, 
and  amounted  to  a  smooth  precession  about  the  Laplace  plane  which  lies 
between  the  perturber's  orbital  plane  and  the  planet's  equatorial  plane. 
The  Laplace  plane  moved  closer  and  closer  to  the  perturber's  plane  as 
the  particle  was  further  away  from  the  planet;  the  plane  was  warped. 

Many  satellite  perturbers  were  not  studied  but  I  would  expect  that  the 
plane  of  symmetry  of  Saturn's  ring  would  be  due  to  the  global  perturba¬ 
tion  of  all  satellites  with  only  small  local  effects  due  to  individual 
sate  Hi  tes  . 

Mico^ni :  It  would  be  interesting  to  see  similar  calculations  done  for 
Venus . 


Cjck:  The  exchange  which  we  have  just  heard  appears  to  neglect  the  aver¬ 
aging  effect  of  integrating  along  the  line-of-sight .  If  we  start  with 
symmetry  about  the  solar  equator  close  to  the  Sun  and  end  with  the 
Gegenscnein  necessarily  being  close  to  the  ecliptic  at  the  other  end, 
we  must  have  local  latitudes  of  symmetry  which  correspond  to  a  steady 
decrease  of  the  apparent  inclination  of  the  plane  of  symmetry  as  we 
move  away  from  the  Sun.  If,  in  addition,  the  local  Laplace  plane  plays 
a  role,  then  the  pattern  will  show  a  wavy  structure  -  an  effect  which 
is  probably  beyond  the  reach  of  observations. 

keir.crt:  Plotting  the  predicted  ecliptic  latitude  of  the  point  of  max¬ 
imum  zodiacal  light  brightness  makes  it  appear  as  if  a  determination  of 
the  plane  of  symmetry  in  the  F-corona  was  not  possible.  This  is  mis¬ 
leading.  The  problem  can  be  avoided  by  plotting  the  inclination  angle 
instead  of  ecliptic  latitude.  Also  it  might  be  useful  to  check  by 
model  calculations  whether  the  approximation  used  is  still  good  at  very 
small  and  very  large  elongations. 

y.ic;oni :  L'se  of  the  inclination  angle  is  a  better  measure  of  the  maximum 
brightness  but  not  of  the  photometric  axis  in  the  F-corona.  We  believe 
that  information  on  the  symmetry  plane  can  best  be  obtained  by  scanning 
perpendicular  to  the  ecliptic  plane  and  by  finding  the  displacement  of 
the  photometric  axis  in  ecliptic  latitude.  Our  method  of  prediction  is 
also  reliable  at  small  and  large  elongations. 
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An  observational  and  theoretical  study  of  the  diffuse 
astronomical  background  sky  brightness  (background  starlight)  is 
carried  out.  The  brightness  is  determined  over  95%  of  the  sky 
using  Pioneer  10  photometric  measurements  in  sky  regions  where  the 
zodiacal  light  is  negligible  (heliocentric  distances  ~  3  A.U.). 

O 

Brightness  levels  are  presented  at  blue  (3950-4850A)  ana  red 

O 

(5900-6800A)  wavelengths.  The  B-R  color  index  distribution  is 
established  over  the  celestial  sphere. 

Pioneer  10  results  are  compared  with  previous  star  count  and 
ground  based  photometric  studies  to  separate  background  starlight 
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into 


its  constituent  parts:  integrated  starlight,  diffuse  galactic 
light  (DGL) ,  and  cosmic  light.  Significant  errors  are  found  in 
published  star  count  results  at  low  galactic  latitudes.  The 
galactic  latitude  (b^)  and  longitude  (1^)  dependences  of  inte¬ 
grated  starlight  and  the  variation  of  DGL  with  b^~  are  determined. 

A  value  of  1.35  t  1.3  3  (V)  is  deduced  for  the  cosmic  light 
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brightness  at  blue  wavelengths  near  the  galactic  poles.  The  inte¬ 
grated  light  iron  discrete  galaxies  adequately  explains  this 
component  of  the  background  starlight. 

A  computer  model  of  our  galaxy,  incorporating  detailed  inter¬ 
stellar  absorption  data  and  consistent  with  the  brightness 
distribution  observed  from  Pioneer  10,  has  been  constructed  for 
the  purpose  of  studying  background  starlight.  The  model  computes 
the  emitted  and  observed  brightness  from  any  volume  element.  Knowl 
edge  of  the  brightness  and  absorption  distributions,  coupled  with 
the  Henvev-Greens tein  phase  function,  permits  calculation  of  DGL 
brightnesses.  The  relation  between  DGL,  direct  starlight,  and 


extinction  along  the  line  of  sight  is  established.  Possible  values 


for  tne  aioeao  and 


for  the  interstellar  dust  at  blue 


wavelengths  are  deduced. 


A  small  influence  of  galactic  spiral  structure  on  surface 
brightness  levels  exists  at  1  =  2°  (Sagittarius  arm),  1  =  69° 

(Orion  arm),  and  1  =  288°  (Carina  arm).  No  indication  of  spiral 
structure  effects  in  the  directions  of  nine  other  brightness  maxima 


and  minima  is  discerned.  A  correlation  between  brightness  features 
and  interstellar  extinction  is  established.  Analysis  of  the  dis¬ 
tribution  of  brightness  peaks  and  dips  and  the  location  of  the 
brightness  symmetry  plane  establishes  the  sun's  elevation  above  the 
galactic  plane  as  12.2  i  2.1  pc. 
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Abstract 


Asano  and  Yamamoto ' s (1975)  solution  of  the  scattering  by  homogeneous 
spheroids  of  size  comparable  to  the  wavelength  has  been  checked  and  used 
to  calculate  selected  results.  The  entire  solution  has  been  rederived  as 
a  check  and  presented  in  detail  as  a  guide  to  further  development.  Test 
results  between  the  computer  code  developped  here  and  that  of  Asano's 
show  good  agreement;  in  addition,  results  from  two  other  recently  developped 
theoretical  methods  agree  with  the  corresponding  results  here  to  within 
published  accuracy.  Comparisons  are  also  shown  to  laboratory  measure¬ 
ments;  those  of  the  Microwave  Scattering  Facility (State  University  of  New 
York — Albany,  now  at  the  University  of  Florida  in  Gainesville)  agree  well 
for  all  measurements.  Finally,  results  are  shown  for  the  scattering  by 
a  selection(cf .  Table  I)  of  randomly  oriented  ensembles  of  spheroids. 

All  six  independent  elements  of  the  Mueller  matrix  along  with  the  extinction, 
scattering,  and  phase  lag  cross-sections  and  the  asymmetry  factor  have 
been  calculated.  From  the  calculations  here  and  previous  results  reported 
by  other  authors,  the  scattering  by  spheroids  has  been  discussed  and  charac¬ 
terized. 
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Abstract 


A  hypothetical  selfconsistent  picture  of  evolution  of  prestellar  interstellar 
dust  through  a  comet  phase  leades  to  predictions  about  the  composition  of  the 
circum-solar  dust  cloud.  Scattering  properties  of  thus  resulting  conglomerates 
with  a  "bird  's-nest'1  type  of  structure  are  investigated  using  a  micro-wave 
analogue  technique.  Approximate  theoretical  methods  of  general  interest  are 
developed  which  compare  favorably  with  the  experimental  results.  The  principal 
features  of  scattering  of  visible  radiation  by  zodiacal  light  particles  are 
reasonably  reproduced.  A  component  which  is  suggestive  of  &-meteoroids  is  also 
predicted . 
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debris  —  interplanetary  dust  —  zodiacal  light. 


